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In the first three years since the discovery of Fe-based high T c superconductors, scanning tunneling 
microscopy (STM) and spectroscopy have shed light on three important questions. First, STM 
has demonstrated the complexity of the pairing symmetry in Fe-based materials. Phase-sensitive 
quasiparticle interference (QPI) imaging and low temperature spectroscopy have shown that the 
pairing OP varies from nodal to nodeless s=b within a single family, FeTei- x Secc. Second, STM has 
imaged C4 — >• C2 symmetry breaking in the electronic states of both parent and superconducting 
materials. As a local probe, STM is in a strong position to understand the interactions between these 
broken symmetry states and superconductivity. Finally, STM has been used to image the vortex 
state, giving insights into the technical problem of vortex pinning, and the fundamental problem 
of the competing states introduced when superconductivity is locally quenched by a magnetic field. 
Here we give a pedagogical introduction to STM and QPI imaging, discuss the specific challenges 
associated with extracting bulk properties from the study of surfaces, and report on progress made 
in understanding Fe-based superconductors using STM techniques. 

PACS numbers: 68.37.Ef, 74.25. Jb 
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I. INTRODUCTION 

The 2008 discovery of high-T c superconductivity in 
Fe-based materials marked an exciting turning point in 
the study of unconventional superconductivity pQ. Until 
then, cuprates were an anomalous island, hosting a zoo 
of confusing properties, with little outside perspective 
to determine which of these properties deserved deeper 
theoretical and experimental attention. Furthermore, 
the cuprates' brittle material properties, gross electronic 
anisotropy, and poor normal state electrical conductivity 
made them challenging to incorporate into widespread 
technology. The Fe-based superconductors provided a 
foil for comparison and a fresh start for all involved in 
the cuprate quagmire. Furthermore, their relative mal- 
leability, isotropy, and metallic normal state led to high 
hopes for useful applications. These considerations have 
launched a whirlwind of research. Three years into this 
new exploration, we review their properties, focusing on 
scanning tunneling microscopy (STM) experiments. 

STM has proven to be an ideal tool to study cor- 
related electron materials. These materials are prone 
to nanoscale inhomogeneities whose effects may broaden 
spectral features or transitions measured by bulk tech- 
niques. STM has been applied with great success to 
cuprate superconductors, addressing pairing symmetry, 
gap inhomogeneity, dopant placement, vortex pinning, 
and competing phases [2]. In recent years, STM has 
also been used to gain insight into momentum space, via 
quasiparticle interference (QPI) imaging. QPI imaging 
can even provide a phase-sensitive determination of the 
superconducting order parameter (OP). 

An early review of STM of Fe-based superconductors 
was written by Yin et al[3 . More recent, comprehen- 
sive reviews of the thousands of papers to date on Fe- 
based superconductors have been written by Johnston [4 
and Stewart [5]. These longer reviews include short sum- 
maries of STM results. Here we give a more thorough 
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FIG. 1. Schematic of a scanning tunneling microscope. A 
voltage is applied between the conducting sample surface and 
a sharp metallic tip, leading to a measurable tunneling current 
whose magnitude decays exponentially with tip-sample sepa- 
ration. In topographic mode, the measured tunneling current 
/meas is compared to the setpoint current J se t, and the differ- 
ence (the error signal) is fed back to the z piezo to control the 
tip height. 




FIG. 2. Schematic of tip-sample tunneling. Energy is along 
the vertical axis, and DOS of the sample and tip are shown 
along the horizontal axes. Filled states are shown in green. 
In this case, a positive bias voltage V has been applied to 
the sample, which effectively lowers its Fermi level by energy 
eV with respect to the Fermi level of the tip. This allows for 
filled states from the tip (right) to tunnel into empty states 
in the sample (left). The tunneling current is measured by an 
external circuit. 



review including a pedagogical introduction to STM, a 
discussion of surface configurations, and several signifi- 
cant new results. 

We start with an introduction to STM in section [TTJ in- 
cluding a pedagogical explanation of QPI. In section [HI 
we give a brief overview of Fe-based superconductivity. 



In section \LV\ we discuss the crystal structure and sur- 
face characterization of the several families of Fe-based 
superconductors. In section [V] we report on measure- 
ments of the superconducting OP, focusing on phase sen- 
sitive measurements via QPI imaging. In section |VT] we 
discuss the parent compound and competing electronic 
orders. In section [VlTJ we discuss the vortex state. In sec- 
tion VIII we conclude with suggestions for future STM 
experiments that could shed additional light on these ma- 
terials. 



II. SCANNING TUNNELING MICROSCOPY 

A scanning tunneling microscope consists of a sharp 
metallic tip which is rastered several Angstroms above an 
electrically conducting sample surface [6]. The position of 
the tip can be varied with sub- A precision by means of 
a piezoelectric scanner. When a voltage V is applied be- 
tween tip and sample, a current will flow. This current 
can be measured as a function of (x, y) location and as 
a function of V. The microscope is illustrated schemat- 
ically in figure [I] There are several excellent textbooks 
about STM[7H9]. 



A. Tunneling Current 

When a positive voltage V is applied to the sample, 
with respect to a grounded tip, the Fermi level of the 
sample is lowered with respect to that of the tip. Elec- 
trons will flow primarily from the filled states of the tip 
into the empty states of the sample, as illustrated in fig- 
ure [2] The elastic tunneling current from the tip to the 
sample, for states of energy e with respect to the Fermi 
level £p = of the sample, is given by 



-2e-^|M| 2 - 



We) 

Pt (s-eV)f(e-eV) ■ Ps (s) [1 - /(e)] 



(1) 



# filled tip states 
for tunneling from 



# empty sample states 
for tunneling to 



where there is a factor of 2 for spin, — e is the electron 
charge, |M| 2 is the matrix element for the tunneling bar- 
rier, p s (e) is the density of states (DOS) of the sample, 
p t (e) is the DOS of the tip, and f(e) is the Fermi distri- 
bution, 
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(2) 



Though the dominant flow of electrons for positive 
sample voltage V will be from tip to sample, there will 
also be a smaller flow of electrons from sample to tip, 
given by 
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We) = -2e-— |M| 2 - 

p s (e)/(e) • Pt (e-eV)[l-f(e-eV)}. (3) 

V V 

# filled sample states # empty tip states 

for tunneling from for tunneling to 

Summing these counter-propagating currents, and inte- 
grating over all energies e, gives a net tunneling current 
from tip to sample 

Aire f°° 

It^s = —^j \M\ 2 p s (e)p t (e-eVy 

[f(e-eV)-f(e)]de. (4) 

This expression can be simplified at low temperature, 
where the Fermi function cuts off very sharply at ep • (For 
example, at T=4.2K, the cutoff width is ~ 2k B T = 0.72 
meV, which can be compared to a typical superconduct- 
ing gap in Fe-based superconductors of ~ 5 meV.) In the 
approximation of a perfectly abrupt cutoff, the integrand 
is negligible except in the range < e < eV. (Likewise, a 
negative bias voltage — \V\ applied to the sample, would 
give an integration range of — \eV\ < £ < 0.) The tun- 
neling current for positive V is therefore 



The matrix element \M\ 2 can be calculated for an 
s- wave tip [11 . The tunneling probability through the 
vacuum barrier is approximated by \M\ 2 = e~ 2j with 
7 = d^2m(p/h 2 , where m is the mass of the electron, d 
is the width of the barrier (tip-sample separation) , and (p 
is the height of the barrier, which is some mixture of the 
work functions of the tip and sample, typically ~ 3 — 5eV 
for metals (see table |l|. In summary, the tunneling cur- 
rent is fairly well approximated by 

Aire h f*^ f eV 
it-, « --jr-e-'y/^-fHiO) J p,(e)de. (8) 



Element 


<P (eV) 


Element 


if (eV) 


Element if (eV) 


Se 


5.9 


W 


4.32-5.22 


As 3.75 


Pt 


5.12-5.93 


Te 


4.95 


Ca 2.87 


Ir 


5.42-5.76 


Fe 


4.67-4.81 


Sr 2.59 


Au 


5.31-5.47 


Bi 


4.34 


Ba 2.52 


Co 


5.0 


Pb 


4.25 


K 2.29 



TABLE I. Work functions of several relevant elements. 
Ranges correspond to different crystalline orientations [T2]. 



Aire f eV 

h^s « J \M\ 2 p s (e)p t (e - eV)de. (5) 

In reality, equation [5] will be modified by an apparent 
smearing of energy features with width ~ Ak B T (com- 
posed of ~ 2k B T for the sample and ~ 2k B T for the 
tip). 

A second simplification comes from the choice of a tip 
material with a featureless DOS near ep. Typical choices 
are W, Pt, or Ptlr. If the tip DOS is fiat in the energy 
range of interest, then p t (e + eV) can be treated as a 
constant and taken outside the integral, giving 

Aire P eV 
It^s ~ — R-Pt(0) ^ \M\ 2 p s (e)de. (6) 

The final simplification is due to Bardeen's demonstra- 
tion that under several realistic assumptions, the matrix 
element for tunneling will be virtually independent of the 
energy difference between the two sides of the barrier [TO]. 
In particular, the matrix element will remain unchanged 
even if one side transitions from the normal state to the 
superconducting state. To a reasonable approximation, 
the matrix element can therefore be taken outside the 
integral, giving 

Aire o f eV 
I^s^-^\M\ 2 p t (0)j^ p s (e)de. (7) 



B. Topography 

In topographic mode (illustrated in figure [I]) , the tip- 
sample bias voltage is fixed, and a feedback loop holds 
the tunneling current constant by varying the tip height 
while it is scanned over the sample surface. Therefore, 
the tip follows a contour of constant integrated DOS 
(DOS). In the case of a homogeneous metal, the contour 
of constant DOS corresponds to the geometric topogra- 
phy of the sample surface. However, if the local DOS 
varies spatially, the resulting image contains a mixture 
of DOS and true topographic information. When the 
tip-sample bias voltage can be set far from the energy 
range of spatially inhomogeneous states, the contribution 
of the geometry dominates the topographic image, as de- 
sired. The recorded tip height is linear in the geometric 
topography, but logarithmic in the integrated DOS. 



C. Spectroscopy 

In addition to revealing the geometrical surface struc- 
ture of a sample, STM can also measure the sample DOS 
as a function of energy, up to several eV from the Fermi 
level in both occupied and unoccupied states (where the 
upper bound is set either by catastrophic surface destruc- 
tion at high fields, or by breakdown of the tunneling ap- 
proximation). This is accomplished by sweeping the bias 
voltage V and measuring the tunneling current I while 
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maintaining constant tip-sample separation d. By differ- 
entiating I(V), the conductance dl/dV is found to be 
proportional to the sample DOS. 




p s {e)de 



(9) 



To reduce noise, a lock-in technique is typically used: a 
small AC modulation is summed with the DC bias volt- 
age, and the resultant tunneling current is demodulated 
to yield dl/dV. The dl/dV spectrum is a good mea- 
sure of the sample DOS up to an overall constant which 
depends on the tip-sample bias voltage and current set- 
points. If these dl/dV spectra are recorded at a dense 
array of locations in real space, spatial variation in the 
sample DOS can be extracted. 



D. QPI 

When materials are inhomogeneous and electrons 
are highly correlated, a full understanding of their 
properties requires knowledge of both the real space and 
momentum space behavior of quasiparticles. QPI (QPI) 
imaging provides some of this complementary knowl- 
edge. When quasiparticles scatter off defects or other 
structures within the crystal, energy-dependent standing 
waves form. (The special case of standing waves at 
the Fermi energy are known as Friedel oscillations [13].) 
The resulting interference patterns in the quasiparticle 
DOS can be imaged with scanning tunneling microscopy 
and spectroscopy [T4HT7] . The Fourier transform of the 
real-space interference patterns highlights the dominant 
sets of quasiparticle moment a [18-23^. This combination 
of imaging and analysis is often called Fourier-transform 
scanning tunneling spectroscopy (FT-STS). FT-STS is 
a powerful technique because it simultaneously yields 
energy-dependent real-space and momentum- space in- 
formation on the quasiparticle wavefunctions, scattering 
processes, and coherence factors. This information can 
be used to distinguish between candidate superconduct- 
ing OPs[24H27]. 

Metals In an ideal metal, the Landau quasiparticle 
eigenstates are Bloch wavefunctions characterized by 
wavevector k and energy e. Their dispersion relation, 
£(fc), can be measured with momentum resolved tech- 
niques such as angle resolved photoemission spectroscopy 
(ARPES)[29 j. By contrast, real space imaging tech- 
niques, such as STM, cannot directly measure e(k). The 
local DOS, LDOS(^,r), is related to the fc-space eigen- 
states ipk(r) by 



LDOS(£, f) oc \Mr)\ 2 S(E - e(k)) (10) 
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ip k (r) = sin(^-^ f) I 2 I 
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FIG. 3. Schematic of QPI. (a-b) BZs for two hypotheti- 
cal materials show that large JDOS can arise in two different 
ways. The thin blue lines represent CCEs. (a) Nested re- 
gions, which will give large JDOS, are emphasized by thicker 
blue lines, (b) Regions of shallow dispersion, where the CCEs 
are far apart, will also give large JDOS. Examples of these 
so-called 'hotspots' are noted with blue dots. In both (a) and 
(b) , pairs of eigenstates k\ and are nested with wavevector 
q — k\ — ki. (c) Schematic showing how two different spatially 
homogeneous Bloch states (left) can mix to form modulations 
of the wavefunction (middle) and the LDOS (right). For 
superconducting Bogoliubov quasiparticles, the filled states 
(bottom right) and empty states (top right) of the LDOS are 
spatially out of phase with each other [28], 



Substitution of a Bloch wavefunction into equation [l0| 
shows that LDOS(£?,r) does not carry any directly ob- 
servable spatial modulation at wavevector k. 

Spatial structures such as impurities, crystal defects, 
magnetic vortices, or static or fluctuating spatial orders 
cause elastic scattering which mixes eigenstates of dif- 
ferent k but the same e(k). A full treatment of this 
mixing was carried out using a Green's function ap- 
proach for a normal metal in the presence of a weak local 
potential [30]. A simpler picture can be explained as fol- 
lows: elastic scattering mixes states that are located on 
the same quasiparticle contour of constant energy (CCE) 
in /c-space. For example, the blue lines in figures |3p and 
b show CCEs for two different hypothetical materials. 
When scattering mixes states k\ and &2, the result is a 
standing wave in the quasiparticle wavefunction ^ of 
wavevector [k\ — k^jl. Since LDOS is proportional to 
l^/cl 2 , the LDOS will contain an interference pattern with 
wavevector q = k\ — &2, as sketched in figure [3J3. LDOS 
modulations of wavelength A = 2ir / 'q can be observed by 
STM as spatial modulations of the differential tunneling 
conductance dl/dV. 

Although the full Green's function treatment includes 
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additional terms, the experimentally observed amplitude 
of these modulations can be simply understood from 
Fermi's golden rule: 

2-7T -> 

w(i ^ f) <x —iViqj^miE^k^nfiEf,^) (11) 

where Ei = Ef for elastic scattering, q = kf — k^ V{q) 
is the Fourier component of the scattering potential at 
wavevector q, and Ui and n f are the density of the initial 
and final states. 

It is apparent from equation [TT] that a particular 
wavevector q can dominate the QPI at energy E, if 
the CCE contains a large joint DOS (JDOS) of /c-pairs 
connected by that q. Large JDOS for a particular q can 
arise from nested regions where two equal-energy con- 
tours are roughly parallel and separated by <f, as shown 
in figure [3^i, or from large flat regions in /c-space, where 
DOS(£) oc 1/| Vfe (E)\ is large, i.e. the CCE are farther 
apart, as shown in figure [3}d. These relatively flat regions 
of /c-space are often called 'hotspots'. (The extreme case 
of shallow dispersion is a van Hove singularity (vHS). 
A vHS has large JDOS with all other fc-states, thus 
it will give rise to a FT-STS image which mimics the 
original band structure [31].) The simple JDOS picture 
can be experimentally verified for a given material by 
a careful comparison between the autocorrelation of 
the ARPES-derived band structure (i.e. the JDOS), 
and the Fourier transform of the STM-derived QPI 
modulations [32H35] . 

Superconductors The Bogoliubov quasiparticles in a 
superconductor are also Bloch states, but with dispersion 

£±(£) = ±\A(fc) 2 + |A(fc)|2 (12) 

where |A(fc)| is the ^-dependent magnitude of the super- 
conducting energy gap at the Fermi surface (FS). (The 
FS is the CCE for e(k) = in the normal state.) Elas- 
tic scattering of Bogoliubov quasiparticles can also result 
in conductance modulations, which have been treated in 
detail by several authors [30 ], l36ti42] . In a superconduc- 
tor there is an additional complication from coherence 
factors [36] l38| [40-42 . In the simple picture, Fermi's 
golden rule is augmented as 

w(i ->• /) a —lukXk^Vk.vl/lV^niiEuk^nfiEf^f) 

(13) 

where the plus sign is for magnetic scatterers, the minus 
sign is for non-magnetic scatterers, and l^l 2 and \vk\ 2 are 
the probabilities that a pair of states with wavevectors 
±fc is empty or filled, respectively [43]. These Bogoliubov 
coefficients Uk and Vk are given by 

A(fc) [l 7 e(k) \ r 




FIG. 4. Determination of the cuprate d-wave OP via QPI 
imaging, (a) Schematic of the cuprate BZ. The normal state 
FS is shown in red. As the material becomes superconduct- 
ing, a nodal gap opens, resulting in the blue contours of con- 
stant energy (CCEs). The DOS is largest in the flat regions 
at the end of the banana-shaped CCEs. The QPI signal is 
dominated by the 7 ^ 's connecting this octet of high-DOS 
regions, (b-c) Occupation factors Uk and Vk for a d-wave gap, 
mapped in the same BZ as in (a), (d) Z(r, E — 4.4 meV) im- 
age of Ca2- cc Na cc Cu02Cl2 at zero field and T=1.6K (setup: 
Sample = -100mV;/ se t = 100 pA). (e) FT-STS image: 8- 
fold-symmetrized Fourier transform of (d), showing the 7 q^s 
expected from the connections between the 8 'hotspots' in 
(a), (f) The difference between high-field and zero-field QPI: 
Z(q, E, B = 11 T) - Z(q, = T) for £=4.4 meV. This 
figure is taken from Ref. [26] 



Therefore, in a superconductor, a large JDOS will result 
in large QPI signal only when the additional constraint 
is satisfied that the coherence factor is non-zero. 

Hanaguri[28] has also pointed out that for supercon- 
ductors, the QPI dispersions qj(E) are particle- hole sym- 
metric, i.e. qj(E) = qj(—E) where qj represents a dom- 
inant scattering wavevector as exemplified in figure [3^i 
and b, and E is measured with respect to £p. However, 
for a given qj , the real-space conductance modulations at 
+E and — E are spatially exactly out of phase with each 
other, as shown in the right panel of figure The QPI 
signal can therefore be enhanced relative to the back- 
ground by computing the ratio map 

The ratio also eliminates possible spurious effects from 
the tunneling matrix element [28 . 

Example: Cuprate Superconductors Although the 
pairing symmetry of cuprate superconductors has long 
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been known [44) [45] . a recent STM study on supercon- 
ducting Cas-^Na^CuC^C^ (x -0.14 and T c -28K) 
demonstrated the use of QPI imaging as a phase-sensitive 
probe of a superconducting OP [26 . Here we will walk 
through this example as preparation to better understand 
the determination of the FeSe OP using QPI imaging dis- 
cussed later in section IV1 

In optimally doped cuprates, there is a single Cu d 
band crossing the Fermi level, resulting in a single large 
hole pocket centered at the M point of the single-CuCV 
plaquette Brillouin zone (BZ). The normal state FS is 
shown schematically as a red line in figure [4^i. Below 
T c , a gap opens with four nodal points, resulting in the 
banana-shaped CCEs shown schematically in blue in fig- 
ure [4^i. The JDOS is therefore dominated by the 7 in- 
equivalent ^-vectors connecting the octet of flat regions 
at the ends of the 'bananas' [2j 136} 146]. The measured 
FT-STS image in figure [4^ indeed shows all 7 of these 
expected dominant ^-'s. 

The dominant superconducting OP in cuprates is d- 
wave, with the gap changing sign along the BZ di- 
agonals, as denoted by the white and gray shading in 
figure [4^i. The sign change in results in sign-changing 
Bogoliubov coefficients Uk and Vk, shown in figure [4}> 
c. The 7 q^s therefore have different coherence factors 
| u* kf ±Vk { v% f | 2 . The question is, within a given sample, 
which type of scattering dominates, and therefore which 
sign applies within the coherence factor? For weak scalar 
potential scattering, the coherence factor is suppressed 
for <£/'s that preserve the sign of namely <?4, and 
q$. For scattering from magnetic impurities, the coher- 
ence factor is suppressed for <^'s that change the sign of 
Afc, namely <?2, <?3, qs-, qi- When a magnetic field is ap- 
plied, the pinned vortices serve as an additional form of 
gap-sign-preserving scattering [41 , thus suppressing qi, 
Q3i Q6, Q7 while preserving or enhancing gi, #4, and q$. 
This effect is observed experimentally in figure [4f. The 
phase of the superconducting gap in cuprates is thus re- 
vealed by the evolution of the QPI upon application of 
the magnetic field. 



III. FE-BASED SUPERCONDUCTOR 
OVERVIEW 

The common feature of all known Fe-based supercon- 
ductors is an FeX square lattice, where X may be a pnic- 
togen (As or P), or a chalcogen (Te, Se, or S). Between 
these FeX layers there may be an additional oxide layer 
('1111' materials), ionic layer ('122' or '111' materials), 
more complex intermediaries (e.g., '42622' materials), or 
nothing ('11' materials). Representative crystal struc- 
tures are shown in figure [5] 

Fe-based superconductors, like cuprates, exhibit a 
dome-shaped phase diagram with an antiferromagnetic 
(AF) parent compound. Additionally, they undergo a 
structural phase transition which is often but not always 
simultaneous with the magnetic phase transition. Phase 
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FIG. 5. Structures of 4 families of Fe-based superconductors, 
(a) 1111 (specifically, LaFeAsOi-.F.), from Ref.gTl (b) 122 
(specifically, BaFe2As2), from Ref . [48l (c) 111 (specifically, 
LiFeAs), from Ref.Hfl (d) 11 (specifically, FeTe), from Ref. EH 



diagrams representative of three different families of Fe- 
based superconductors are shown in figure [6| 
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FIG. 6. Phase diagrams of several Fe-based super- 

conductors, representative of the 1111, 122, and 11 fam- 
ilies, (a) SmFeAsOi-.F. [51]; (b) LaFeAsOi-.F. [52]; 
(c) Ba(Fei_ a ,Co a ,) 2 As 2 [53]; (d) FeTei_.Se. 54 . The or- 
thorhombic transition marked in (d) comes from x-ray diffrac- 
tion measurements on a powder sample with stoichiometry 
Fei.oiSe [55]. In each diagram, stars mark the phase space 
locations of STM experiments to be discussed later. 



There are a number of points of confusion when con- 
sidering the real space unit cell or the BZ (BZ) of these 
materials. The common FeX layer in all of these materi- 
als (X=As,P,Se,Te,S) should be more accurately written 
as Fe2X2, because half of the X lie above the Fe plane, 
while half lie below, as can be seen in figure [5] In the 
high temperature state, the material is tetragonal, with 
a unit cell twice as large (dashed square in figure [7^) and 
BZ half as large (figure [Tp) as might be naively expected 
from the Fe square sublattice (solid square in figure [7^, 
and BZ in figure up). In the underdoped region of the 
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phase diagram, the material undergoes a tetragonal to 
orthorhombic phase transition as the temperature is low- 
ered. Then the unit cell doubles again (figure [7)3), and 
the BZ halves again (figure^). However, the difference 
between the new ao and bo axes in the orthorhombic 
state is typically < 1% (e.g. see tables A8 and A9 in 
Ref. |4|) which is below the calibration resolution of the 
piezo scant ubes typically used by STM. Therefore the 
ao and bo axes cannot be distinguished by STM unless 
a twin boundary between orthorhombic domains can be 
found [561 [57]. 

There is an additional source of confusion regarding 
the 122 materials. The 1111, 111, and 11 materials have 
only a single FeX layer per unit cell. However, the 122 
materials have a double layer along the c-axis, so their 
unit cell is technically body-centered tetragonal, with a 
BZ that is rotated by 45° from the tetragonal BZ of the 
other Fe-based families. In this review, since STM im- 
ages only a single layer at a time, we will ignore this 
complication and we will talk about the 122 unit cell and 
BZ in the same tetragonal/orthorhombic notation as is 
used for the other Fe-based superconductors. 

Fe-based superconductors are more complicated than 
their cuprate cousins because all five Fe d-bands cross 
the FS (FS), in contrast to the single Cu d band crossing 
the FS in cuprates. The schematic bands are shown in 
the unfolded BZ in Fig [7J3. Local density approximation 
(LDA) band structure calculation showed that Fe-based 
materials are semi-metals with hole bands at the V point 
and electron bands at the M point, in the tetragonal BZ 
(figure [7]i) [58] . The materials may be electron-doped or 
hole-doped, but over a large range of doping there is still 
usually significant nesting between the hole and electron 
FSs. 
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FIG. 7. Real space and momentum space unit cell of Fe- 
based superconductors, (a) Fe2X 2 lattice, where the red dots 
represent Fe atoms, and the blue filled (open) circles represent 
X=As, P, Se, Te, or S atoms above (below) the Fe plane. The 
solid square shows a unit cell with only one Fe; because of the 
two inequivalet X sites above and below the Fe plane, this 
is not the true unit cell of the crystal. The dashed square 
shows the unit cell of the tetragonal structure typically seen 
at high temperatures, (b) The a and b axes of the low T 
orthorhombic unit cell (black rectangle) typically differ by 
<1%; the difference is exaggerated in this figure for visual 
clarity, (c) Unfolded BZ of the unit cell containing only a 
single Fe atom. Band structure calculations typically predict 
as many as three inequivalent hole FSs (blue), centered at 
the r point, and 2 inequivalent electron FSs (red), centered 
at the M points. Variations in structure and charge doping 
can reduce the number of hole FSs, in some cases leaving only 
a single a surface at the T point, (d) Calculated LDA FSs 
in the tetragonal BZ of Ba(Feo.gCoo.i)2As2[59]. Note that 
in this doubled real space unit cell and halved BZ, all hole 
FSs from (c) have now been folded into the central T point, 
while the two electron FSs have been folded onto the same 
M point, (e) Calculated LDA FS in the orthorhombic BZ for 
CaFe2As2 in the AF state [60]. The real space unit cell has 
doubled again, and the BZ has halved again with respect to 
(d). 



IV. SURFACE CONSIDERATIONS 

Like the cuprates, the Fe-based superconductors have a 
layered structure. Consequently, high-quality single crys- 
tals can be mechanically cleaved to obtain atomically flat 
and clean surfaces suitable for characterization by surface 
sensitive probes such as STM (STM) and spectroscopy 
(STS). With effort, samples as small as ~ 100 jam can be 
cleaved and imaged (although typically 1mm is a more 
convenient size). Just because a surface is clean and 
atomically flat does not necessarily mean it is a window 
into the bulk properties of the material. Surfaces may 
have different carrier concentration than the bulk, and 
may exhibit electronic or structural reconstructions due 
to the reduced coordination number. In this section, we 
will discuss the surface characteristics for the different 
families of Fe-based superconductors. 



A. 1111 Materials 

Although LaFePO was found to superconduct with 
T c ~5K in 2006 [61 , the Fe-based superconductors first 
rocketed to worldwide attention in February 2008 with 
the discovery of superconductivity up to T C =26K in 
LaFeAsOi-^F^l . These so-called '1111' materials have 
in common iron-pnictogen planes, separated by and ion- 
ically bound to oxide planes. Because of the relatively 
high T c , and the immediate superficial similarities be- 
tween LaFeAsOi-^F^ and the cuprates (e.g. the lay- 
ered structure, dome-shaped phase diagram, and sur- 
prising presence of Fe suggesting a possible spin-based 
superconducting mechanism), the reaction was immedi- 
ate. Research labs around the world swiftly turned their 
attention to the new materials. Within a few months, 
RFeAsO had been formulated with almost a dozen rare 
earth elements (R) [62-66 , raising T c up to 56K[67j. 
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Large single crystals proved challenging to grow, so 
studies of intrinsic anisotropy or inhomogeneity remained 
out of reach. A few early, painstaking experiments man- 
aged to isolate single crystals as large as ~ 100 jam for 
transport [68J or ARPES [69 . However, early STM ex- 
periments on polycrystalline samples [70], ITT] likely suf- 
fered from degraded surfaces. The 'as-grown' surface of 
a ~100 /im single crystal of SmFeAsOi-^F^ (with nom- 
inal £=0.2, and measured T C =45K, denoted by a red ★ 
in figure]^) could not be imaged with atomic resolution, 
but did show a reproducible ^-shaped gap of the ap- 
proximate expected bulk value, that disappeared around 
the bulk superconducting T c [72]. SmFeAsOi-^F^ spec- 
tra are shown in figures [15^ and|2Tfc, and will be discussed 
further in section IVBI 

In late 2009, a breakthrough allowed the growth of 
large single crystals of the '1111' materials using NaAs 
flux at ambient pressure [73]. Their growth method has 
since been replicated by other groups [74]. The first 
STM study on a cleaved single crystal '1111' material 
LaFeAsOi-^F^ (both parent compound, and nominal 
composition £=0.1, denoted by red Vs in figure J6Jd) ap- 
peared a year later [75]. Zhou et al were able to cleave and 
image atomically flat FeAs and LaO planes (figures 19 
a and b), neither of which showed any structural recon- 
struction of the surface. 

Although '1111' single crystals may cleave beautifully, 
surface-sensitive studies are expected to be problematic 
due to the ionic nature of the interlayer bonding, which 
is expected to result in a polar catastrophe at the cleaved 
surface [76]. The surface electronic structure is likely to 
be different from that of the bulk. In the best case, there 
is just a rigid band shift which affects measurements of 
carrier density and Fermi level [77]. In addition, atomic 
reconstruction, adsorption of charged contaminants, and 
electronic reconstruction might occur at the surface in 
order to compensate for the charge imbalance [78]. In 
fact, detailed calculations predict that several new sur- 
face bands appear in cleaved LaFeAsO i-^F^, on both 
La-terminated and As-terminated surfaces [79]. Zhou et 
al found strong evidence for an electronic surface state 
on the LaO plane: the 2-dimensional electronic standing 
waves were so pronounced that the atomic lattice was 
almost completely concealed (see figure [lffi)) ■ 

The existence of dramatically different surface states in 
1111 materials is further supported by ARPES measure- 
ments of the enclosed FS area on the cleaved surfaces 
of LaFePO [80] and LaFeAsO [ST]. In both cases, the 
ARPES-measured FS area falls short of the bulk elec- 
tron count that would be expected from Luttinger's the- 
orem. An even more dramatic discrepancy is found in 
another ARPES experiment which sees the surface of 
NdFeAsOo.gFo.i to be hole-doped, although the bulk is 
known to be electron-doped [69]. A follow-up ARPES [82] 
study found an additional large hole-like, T-centered FS 
on LaFeAsO, which was determined to be a surface state 
by its lack of k z dispersion, and by comparison with the- 
ory. When the LaFeAsO was electron-doped through 



the introduction of F, the 2-dimensional surface FS was 
seen to develop a large superconducting gap, likely via 
the proximity effect from the bulk superconducting state. 
The surface state superconducting gap may have different 
magnitude and additional symmetry components, which 
may confuse surface investigations of the bulk supercon- 
ductivity. 

Although large single crystals of LaFePO have been 
available since 2008, there have been no STM studies to 
date, likely because of its low T c , and the surface chal- 
lenges associated with all 1111 materials. 



B. 122 Materials 

In June 2008, a second family of Fe-based supercon- 
ductors was discovered, including the same FeAs plane, 
but only a single layer of intervening A (alkaline earth 
metal) ions separating the FeAs layers, giving formula 
AFe2As2 [83 . Particular members of this '122' family 
are referred to as Bal22, Srl22, Cal22, etc. 

The 122 parent compounds become superconducting 
upon the introduction of hole or electron dopants, or 
the application of chemical or physical pressure. For ex- 
ample, BaFe2As2 becomes a hole-doped superconductor 
upon replacement of Ba 2+ by K + [83], or an electron- 
doped superconductor upon replacement of Fe 2+ by Co 3+ 
or Ni 4+ [84] [85], or a pressure-induced superconductor 
upon replacement of As by the isovalent but smaller atom 
P [86], or application of physical pressure [87 . The emer- 
gence of superconductivity through chemical substitution 
directly into the superconducting layer [84H86] is in stark 
contrast with the cuprates, where the substitution of im- 
purity elements for even a small percentage of the Cu 
atoms can destroy superconductivity [88] . 

Because the 122 materials were successfully grown as 
large single crystals soon after their discovery [89 , they 
were the first target of serious STM study among Fe- 
based superconductors. Many papers have reported on 
their surface properties, yet the structure of the cleaved 
surface remains controversial. 

The 122 materials are expected to cleave with FeAs 
layers intact, but A ions may end up on either of the 
two cleaved FeAs surfaces. If the A ions do not divide 
evenly, each resulting surface will be polar, just as in the 
1111 family. If the A ions do divide evenly, the resulting 
half-A surfaces will be nonpolar, and possibly ordered 
into one of several superstructures. Atomically resolved 
images of cleaved 122 surfaces typically show 1/2 of the 
atoms which would be expected for a complete As or 
Ba/Sr/Ca layer. The STM groups who have achieved 
these images are sharply divided into two camps: those 
who believe that the observed structures result from a 
1/2-layer of Ba/Sr/Ca, vs. those who believe that the 
observed structures result from a reconstruction of a 
complete As layer, i.e. all As atoms are present but 1/2 
of these atoms are 'invisible' to STM. In the following 
sections, we will present the arguments of both camps. 
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Ba/Sr/Ca Surface The first atomically resolved STM 
images of a Fe-based superconductor were reported by 
Boyer et al on hole-doped Sri_ x K :E Fe2As2 (T C =32K), 
cleaved at T = 10K and imaged at T = 5.3K[1)]. The 
low temperature cleave exposed a flat surface with aver- 
age atomic spacing ^4A as would be expected for either 
a complete Sr/K or As layer. These images also showed a 
2x1 stripe reconstruction at 45° to the orthorhombic lat- 
tice. Patches of missing atoms showed windows through 
to an underlying square lattice, also with atomic spacing 
~4A. Fi gure ySa shows that the underlying square lattice 
is laterally shifted by half a unit cell from the topmost 
2x1 surface, and lies beneath it by 2.8(4)A, close to the 
expected 2A vertical distance between Sr and As layers 
in the bulk. The observed atomic separation of ~4A in 
both striped region and square lattice region are not con- 
sistent with the Fe-Fe distance. Boyer concluded that the 
upper striped region was a reconstruction of the complete 
Sr/K layer, while the underlying square lattice patch was 
the bare As layer. 

Yin et al reported atomically resolved images of 
optimally electron-doped Ba(Fei_ ;E Co :E )2As2 (nominal 
£=0.1, and measured T c = 25. 3K) which showed a 2 x 1 
reconstruction [9T] (see figure^). They found no step 
edges, but occasional 1/2-period lateral shifts. Upon in- 
verse transforming the four brightest spots in the Fourier 
transform of this topography, a pattern emerged which 
suggested that the surface consisted of half of the Ba 
atoms, arranged into stripes of single atom width, allow- 
ing a glimpse between stripes to the As layer beneath 
[3]. A Ba half-layer would imply a nonpolar surface of 
bulk-like carrier concentration, which is consistent with 
the observation of a ubiquitous superconducting gap, and 
other expected features of bulk superconductivity such as 
magnetic vortices |5T]. 

Hsieh et al reported atomically resolved images of 
the parent SrFe2As2, cleaved at room temperature and 
cooled to T = 40K[92 . The surface (figure^) was only 
partially ordered into a similar 2x1 reconstruction as 
shown by Boyer and Yin. Upon heating to 200K and 
reimaging, the long range order of the 2x1 reconstruc- 
tion was lost. Upon recooling, the surface sometimes 
reordered, but sometimes remained disordered [93]. This 
observation suggests that the ordered 2x1 surface re- 
construction may be a metastable ordering of a mobile 
half-Sr layer, resulting from the low temperature cleaving 
process, easily destroyed upon heating, and not always 
reformed upon recooling. It is harder to imagine that a 
clean, complete As layer would disorder so extremely on 
heating to 200K, and fail to reorder upon recooling. 

Massee et al showed a large number of topographies 
of Ba(Feo.93Coo.o7)2As2, cleaved either at room temper- 
ature or at T < 80K [95] (see figure [8j e and f). At 
low T, they often observed the same ~8 A stripes seen 
previously [90, 91 . Like Yin, they observed occasional 
half-period stripe shifts, resulting in a 'ribcage' structure 
(figure [9]p). They also saw less regular, larger rodlike 
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FIG. 8. Various topographic images of cleaved 122 surfaces, 
presumed to be Ba/Sr-terminated. (a-c) 5x5 nm 2 images of 
cold-cleaved surfaces, (a) Sri_ cc K ;E Fe2As2 (T C =32K) cleaved 
at ^10K and imaged at 5.3K. This is believed to be a nearly- 
complete Sr layer, with a glimpse through to the As layer 
below in one region[90 (setup: V^ampie = — 100mV;/ se t = 
200 pA). (b) Ba(Fe .9Co .i)2As2 (T C =25K) cleaved at -25K 
and imaged at ^6K 91 This surface is believed to show every 
other row of Ba, with the intervening As just barely visi- 
ble beneath[3] (setup: V^ampie = — 20mV;/ S et = 40 pA). (c) 
Ba(Feo.9i5Coo.o85)2As2 cleaved at 120K, imaged at 5K with 
very small junction resistance to allow atomic resolution |94| 
(setup: Sample = 20 mV; J Bet = 2nA). (d-f) 20 x 20 nm 2 im- 
ages of warm-cleaved surfaces, showing greater disorder, (d) 
SrFe2As2, cleaved at room temperature, imaged at 40K[92]. 
(e) Ba(Feo.93Coo.o7)2As2 cleaved at room temperature, im- 
aged at 4.2K[95] (junction resistance Rj ~ 0.75 GO), (f) 
Ba(Feo.93Coo.o7)2As2 cleaved at room temperature, imaged 
at 4.2K[25] {Rj ~ 0.75 GQ). 



features of up to 20A width, which evolved smoothly 
without step edges from the more commonly observed 
8A stripes. In room temperature cleaves, they saw dis- 
ordered 8 A stripes (figure^), similar to those seen by 
Hsieh et al. They also sometimes saw a 5.5A V2x V2 
structure, with meandering antiphase boundaries (fig- 
ure [8]p). In total, Massee observed at least 8 different 
surface structures [95], although 2x1 and y/2 x y/2 were 
found to be dominant [96 . These 8 surface structures are 
sketched in figure [9] It is hard to imagine that a single 
clean As layer would display so many different surfaces 
structures. 

Zhang et al imaged Ba(Fei_ a; Co :c )2As2 with six differ- 
ent nominal dopings: x=0, 0.04, 0.085, 0.10, 0.125, and 
0.16. The latter four were found to be superconducting 
with T C =9K, 25K, 22K, and 9K, respectively [94]. Sam- 
ples were imaged at 5K after cleaving at 120K (above the 
100K at which the extra surface superstructure diffrac- 
tion spots started to disappear, as discussed in the follow- 
ing section, but below the 200K at which the surface long 
range order has fully disappeared [92j [95]). Zhang found 
that a y/2 x y/2 superstructure dominated at all dopings 
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FIG. 9. Sketch of various possible partial-Ba surface 
configurations [95 . Large (small) circles indicate the presence 
(absence) of a Ba atom. 



(figure [8^). Theoretical simulations with a half-Ba ter- 
minating layer showed a strong feature at -200 meV, and 
a gross energy asymmetry with larger empty-state DOS, 
in good agreement with the spectra measured on all six 
samples (see figure |20^i) . Zhang therefore concluded that 
this \/2 x y/2 structure was a half-Ba surface. Note that 
the six dopings span the bulk orthorhombic to tetrago- 
nal transition (see figure [6^), so their ubiquitous \/2 x y/2 
structure cannot be due solely to the bulk orthorhombic 
transition. In only two samples out of many studied, they 
saw patches of the 2x1 stripes. 

The argument for the partial-Ba-terminated surface 
was further supported by measurements of the work 
function, ^[97]- Work functions for relevant elements are 
shown in table [TJ taken from Ref. [12j In dozens of Bal22 
samples studied by Massee et al, representing y/2 x y/2, 
2x1, and other surface organizations, the work function 
was never found to be greater than 1.8 eV, and typically 
found to be ~1.5 eV. On Cal22, Massee also found a 
slightly higher work function, ~1.9 eV, suggesting that 
the surface termination was not As, which would be 
expected to show the same, much higher work function 
in both Bal22 and Cal22. 

As Surface Nascimento et al presented STM images of 
BaFe2As2, cleaved at T=80K, revealing a 5.6A square 
lattice [98 . This is the correct spacing for a y/2 x \/2 or- 
ganization of a half Ba layer. In contrast, the complete 
As layer would form a 3.8A, nearly square lattice. How- 
ever, in combination with low energy electron diffraction 
(LEED) data (next section), Nascimento et al hypothe- 
sized that half of the As were invisible due to a structural 
reconstruction. They pointed out that the surface is 1-6% 
orthogonal (in comparison to the 0.7% orthogonality of 
the bulk), and that this surface orthorhombicity persists 
across the superconducting dome [99]. 

In a second study by the same group [100 , more im- 
ages were presented of the parent Bal22 y/2 x y/2 sur- 
face, sporting some white blobs which could be moved 
with the STM tip, which were identified as sparse re- 



maining Ba (figure [l0^) . Li et al hypothesized here that 
half of the As atoms were invisible due to different spin 
environments. Furthermore, they imaged domain walls, 
as seen zigzagging across the bottom half of figure [To| i. 
Larger images showed that these domain walls always 
formed closed loops. On either side of a domain wall, 
the orthorhombic a and b axes maintained their orienta- 
tion, indicating no orthorhombic twin boundaries. But 
the visible and invisible atoms did switch, suggesting an 
antiphase boundary of the the spin state. These domain 
walls were similar to those seen by Massee, shown in fig- 
ure [8]f. (The prevalence of such antiphase boundaries 
may explain the discrepancy between the large Fe spin 
moments computed from density functional theory, and 
the small spin moments measured experimentally [ 10T] .) 




FIG. 10. Topographic images of cleaved parent 122 sur- 
faces, presumed to be As-terminated. Images are scaled 
to the same unit cell size, (a) BaFe2As2 cleaved at 80K, 
showing a \[2 x \[2 surface reconstruction with a domain 
boundary zigzagging across the bottom of the image [TOO] 
(setup: Sample = 23 mV; J S et = 200 pA). (b) SrFe2As2 show- 
ing 2x1 reconstruction on the left, merging smoothly into 
y/2 x \[2 reconstruction on the right, without step edges 102 
[Rj = 0.8 GQ). The left inset at 4x magnification shows that 
at very low junction resistance (lower left strip, Rj = 10 MQ), 
the stripe appears to be two atoms wide. The right inset at 4x 
magnification shows that half of the presumed As atoms are 
invisible in the \[2 x y/2 reconstruction region (Rj = 1.7 GQ). 
(c) Non-superconducting Ca(Feo.97Coo.o3)2As2 showing a 2 x 
1 surface reconstruction, with a domain boundary in the 
lower left corner (setup: Sample = — 50mV;/ S et = 10 pA), 
and 4x magnified inset in the upper right corner (setup: 
Sample = -5 mV; J set = 100pA)[56]. 

Niestemski et al studied SrFe2As2, cleaved both at 
room temperature and at T = 77K |102j . Like Massee, 
they observed a y/2 x y/2 structure in some areas, which 
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merged seamlessly into a 2 x 1 structure in other areas, 
with no intervening step edges (figure [lOfr). They con- 
cluded that this was a complete As top layer, with half 
of the As atoms invisible under standard imaging con- 
ditions (Rj ~ 0.8 GO.). They gave two arguments, as 
follows. First, in the 2x1 regions, very low junction 
resistance imaging allowed visualization of the expected 
number of atoms in a complete As or Sr layer (see bot- 
tom left edge of first inset in figure (ToJd). Second, unlike 
Boyer's previous work on Srl22[90 , after an extensive 
search over large areas of the surface, Niestemski never 
found a step down to an As layer below, which would 
be occasionally expected if the terminating surface were 
Sr. In fact, Niestemski found some steps up to additional 
rod-like structures, which were suggested to be the rem- 
nants of Sr on top of the As surface. Therefore, Niestem- 
ski concluded that both the y/2 x y/2 and 2x1 observed 
structures were different reconstructions of a complete 
terminating As layer. Both Nascimento and Niestemski 
reported identical dl/dV at visible and 'invisible' atomic 
sites. 

Pan et al used STM to study a series of 
Ba(Fei_ ;E Co ;E )2As2, and reported that the undoped 
samples showed the y/2 x y/2 structure, but on increas- 
ing x, the 2x1 structure became dominant [103 . 

Other Experimental Evidence Other surface studies 
of 122 materials have been carried out using photoemis- 
sion spectroscopy (PES), ARPES, and LEED. 

Hsieh et al presented ARPES data on Sri 22, cleaved at 
10K, which showed a strong feature at the BZ X point [92 
(see figure [7]i for BZ reference). This feature was not 
present in the LDA-calculated band structure, but would 
result from the band folding due to a 2 x 1 surface struc- 
ture. When the sample was heated to 200K and cooled 
to 10K for remeasurement, the band- folding artifact at X 
disappeared, indicating that the 2x1 order was totally 
lost by temperature cycling. 

Van Heumen et al also used ARPES to directly show 
a surface state around -200meV near the V point in 
cold-cleaved Bal22, which disappeared after thermal 
cycling [104 . This -200 meV feature was also seen in the 
STM measurements by Zhang in 6 different samples of 
Ba(Fei_ ;E Co; C )2As2 from x=0 to 0.16[94] (figure [20}i) , by 
YinHljg (figure [20)o), by Niestemski in parent Sr 122 [102] 
(figure |20|3) , and appeared in Zhang's theoretical calcu- 
lations for a half-Ba-terminated surface. 

De Jong et al reported photoemission on the undoped 
parent compound BaFe2As2 (with 7% atomic weight 
Sn impurities), cleaved and measured at room temper- 
ature. Variations in photon energy affect the depth 
probed: they used hv = lAOeV photons to probe a few 
Angstroms, and hv = 3keV photons to probe 10s of nm 
from the surface. With low energy (shallow) photons, 
the binding energy peaks for both the Ba Ad and As 3d 
core states showed additional shifted shoulders due to 
surface contributions. These shoulders were much more 
pronounced for the Ba Ad states than for the As 3d states, 



suggesting that Ba sat at the surface. Furthermore, these 
Ba shoulders were much broader than the As shoulders, 
suggesting that the surface Ba atoms were more disor- 
dered than the near-surface As atoms. 

De Jong et al also compared their Ba analysis to anal- 
ogous work on the Ba Ad core states in the cuprate su- 
perconductor YBa2Cu307_ :E . The surface carrier con- 
centration of cleaved YBa2Cu307_ ;E is known to deviate 
from the bulk, resulting in shift to higher binding energy 
for both cations (Ba 2+ ) and anions (O 2- ). In contrast, 
for Bal22 the shift of the surface states was to higher 
binding energy for cations (Ba 2+ ) and to lower binding 
energy for anions (As 1- ), suggesting that surface dop- 
ing was not the leading cause of the surface state. The 
opposite binding energy shifts for Ba and As in Bal22 
are more reminiscent of GaAs, which is known to have a 
Madelung potential shift at the surface. A Madelung en- 
ergy shift is likely to affect the localized, ionic electronic 
levels, rather than the itinerant near-£^ states. There- 
fore, de Jong concluded that the near-£^ states at the 
surface of Bal22 likely did not differ significantly from 
the bulk. 

Massee et al performed LEED on cold-cleaved Co- 
Bal22 surfaces at T = 17K and observed extra spots 
corresponding to both the y/2 x y/2 and 2x1 reconstruc- 
tions (see schematic in figure [TTj a-d). When the sam- 
ple was heated, the y/2 x y/2 and 2x1 spots started to 
lose intensity at T ~ 100K, and by T ~ 200K they were 
gone (only tetragonal spots remained). But when Massee 
recooled to 17K, none of the superstructure spots reap- 
peared. This was consistent with both Nascimento [98J, 
who cleaved warm and never saw the y/2 x \/2 and 2 x 1 or- 
thorhombic spots to begin with, and with Hsieh [92 , who 
also saw the low T superstructure features in both STM 
and ARPES disappear upon warming. Both Massee and 
Hsieh concluded that cleaving usually leaves some frac- 
tion of Ba layer. If the material is cleaved while cold, then 
Ba may be stuck in any number of different metastable 
arrangements. But if the material is cleaved while warm, 
or warmed up after cleaving, then the Ba revert to the 
lowest energy y/2 x configuration, but without long 
range order. Massee speculated that invisibility of this 
^/2x^/2 superstructure in room temperature LEED mea- 
surements may be explained by the proliferation of an- 
tiphase domain walls exemplified in figure |8f. 

Additional information can be gained by varying the 
energy of the incident electrons, and measuring the re- 
sultant variation in LEED spot intensity, in a technique 
called /y-LEED[T06l[T07]. By comparison to models of 
candidate structures, the /V-LEED data can be used to 
determine the terminating element. The Pendry factor 
is a measure of the reliability of a structural fit; it is con- 
sidered good for Rp ~ 0.2, mediocre for Rp ~ 0.3 and 



bad for R P > 0.5 [108] . 

In contrast to Massee's LEED conclusions, an IV- 
LEED experiment performed by Nascimento et al on 
warm-cleaved Bal22 led them to conclude that the sur- 
face termination is As[98 . In Nascimento's experiment, 
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the LEED pattern reflected only the 3.8A tetragonal unit 
cell (figure [lip), with no hint of the 5.6A a/2 x a/2 struc- 
ture observed by STM in the same material (figure 10 1). 
From the orthorhombic (2,0) spot (tetragonal (1,1) spot), 
Nascimento calculated the Pendry factors for Ba, Fe, and 
As terminations, arriving at Rp= 0.57, 0.45, and 0.24, 
respectively. 

One caveat to Nasicmento's work is that the samples 
for LEED study were cleaved at room temperature, be- 
fore cooling to 20K, so the structure from the low-T- 
cleaved surfaces studied by STM in the same paper was 
possibly very different from the room-T-cleaved surfaces 
studied by LEED. This may explain the absence of the 
expected orthorhombic (1,0) spot in the LEED pattern. 

Van Heumen et al investigated the cold-cleaved sur- 
face of BatFei^Co^Ass (£=0.05 and 0.085) via IV- 
LEED [T04] . The LEED images showed fractional spots 
corresponding to both a/2 x a/2 and 2x1 surface struc- 
tures (figure pT]p), which can be explained by the coexis- 
tence of both types of domains. Furthermore, the energy 
dependence of each of these spots was studied from 100 to 
400 eV, and compared with theoretical simulations. One 
set of simulations assumed a 1/2 Ba terminating layer, 
and allowed the top four layers to relax (^Ba-As-Fe2- 
As), resulting in the low Pendry R factors of R p = 0.19 
for a/2 x a/2 and R p = 0.29 for 2 x 1. A second set of 
simulations assumed a terminating As surface, also al- 
lowing the top three layers to relax, but this resulted in 
R p = 0.42 for a/2 x a/2 and R p = 0.48 for 2x1. No- 
tably, in contrast to Nascimento's work, van Heumen's 
work compared calculated IV curves only to the frac- 
tional spots, i.e. those spots corresponding directly and 
exclusively to the a/2 x a/2 or 2 x 1 surface structures 
which were observed in most images by both surface ter- 
mination camps. 

Combining STM and these other measurements 
leads to a picture in which cleaving leaves ~l/2 of the 
Ba/Sr/Ca atoms atop a reconstructed As layer. Possible 
arrangements of the terminating Ba/Sr/Ca atoms are 
sketched in figure (9T 



(a) Real space 



Theory The experimentalists do not agree on the sur- 
face termination, so what do the theorists say? Gao et al 
presented electronic structure calculations on the (001) 
surface of three A122 compounds (A=Ba, Sr, or Ca)[109 . 
They found that it is most energetically favorable for A 
to divide evenly between the two cleaved surfaces, and for 
the half-A terminating layers to organize into superstruc- 
tures as shown in table |TTJ They found that the electronic 
states imaged on surfaces with the a/2 x a/2 structure 
should be representative of the bulk states, while there 
may be some additional surface states near the Fermi en- 
ergy in the 2x1 superstructure. Gao commented specifi- 
cally on the apparent dimerization of the stripes observed 
in Srl22 and given by Niestemski as a primary piece of 
evidence for a complete As terminating surface[102 . Gao 
instead ascribed this apparent dimerization to hybridiza- 
tion between a terminating half-Sr layer and a complete 
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FIG. 11. (a) Real space atomic configuration and (b) ex- 
pected q space LEED pattern, for a y/2 x y/2 arrangement 
of a 1/2-Ba layer on top of a complete As-Fe2-As layer (the 
lower As atoms are not shown), (c) Real space atomic config- 
uration and (d) expected q space LEED pattern, for a 2 x 1 
arrangement of a 1/2-Ba layer on top of a complete As-Fe2-As 
layer (the lower As atoms are not shown). In both (b) and 
(d), the LEED spots are labeled in the tetragonal notation 
(red) and the orthorhombic notation (blue) . (e) LEED image 
of the surface of BaFe2As2, cleaved at room temperature and 
measured at T=20K|98]. Blue circle shows the missing (1,0) 
orthorhombic spot, which would be expected for the y/2 x y/2 
surface observed by STM in the same paper, (f) LEED image 
of the surface of Ba(Feo. 95000.05)2 AS2 (in the bulk orthorhom- 
bic phase), cleaved and measured at T <20K [T0l] . Red circles 
show diffraction spots corresponding to y/2 x y/2, 2x1 and 
1x2 orders which coexist within the same macroscopic elec- 
tron beam spot size. LEED image of Ba(Feo.9i5Coo.o85)2As2 
in the bulk tetragonal phase (not shown here) similarly shows 
all 2 x 1 and y/2 x y/2 spots [104]. 



As layer beneath. 





Bulk structural phase 




AF orthorhombic 


tetragonal 


Bal22 
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TABLE II. Energetically favorable arrangement of the 1/2- 
A terminating surface layer, as calculated by the plane- wave 
basis method 109 . 



Gao suggested that a metastable As termination may 
result from a fast cleave, so he investigated scenarios 
which might lead to half of the As atoms being invisible, 
as claimed in the y/2 x y/2 surface structure of parent 
Bal22 



[98J [100] and Srl22 [102]. Gao found that in 
the Bal22 orthorhombic phase, a small buckling of the 
As terminating surface does strongly affect the STM 
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images, leading to half of the As being invisible. How- 
ever, he found that no such reconstruction occurs in the 
tetragonal phase of Bal22, or in either phase of Srl22. 
Regardless of the presence or absence of As-buckling, 
Gao expected that the low energy electronic structure of 
any As-terminated surface would be strongly modified 
from the bulk. 

Wrap-up To round out this story of STM-measured sur- 
face structure, it is worth mentioning that STM images 
of very lightly Co-doped (non-superconducting) Cal22 
also displayed a 2 x 1 structure (figure [To]:) [56] . match- 
ing Gao's prediction[109 . In the Cal22 experimental 
work, the authors identified the terminating layer as As, 
based on solely on evidence by Nascimento and Niestem- 
ski. However, given newer LEED data[104 , the authors 
have revised their opinion and believe that their termi- 
nating surface is more likely Ca |110| . It may also be 
worth noting that the unambiguously identified cleaved 
FeAs surface of LaFeAsOi-^F^ shows the ~4A lattice 
of the full As layer with no reconstruction and no miss- 
ing atoms [75] (figure 19 1). Neither does the analogous 
Te/Se-terminated surface of FeTei-^Se^ show any recon- 
struction or missing atoms [27j [57l [95j 1111] (figure [T^). 

Table |III| displays a summary of the many studies of 
122 surfaces. Even if we restrict ourselves to the inter- 
pretations of one of the two competing camps, we must 
conclude from the raw data (figures [8] and 10) that there 
are many possible terminating surfaces, close in energy, 
which can result from very slight differences in cleaving 
and temperature history. This conclusion, admittedly 
unsatisfying in its complexity, is supported by Gao's the- 
ory. Despite strong claims by both the half-A termina- 
tion supporters and the As-termination supporters, we 
must ask: it worth debating about these surfaces at all? 
The only relevant question is: which, if any, of these sur- 
faces are representative of the bulk electronic structure 
of the material? 



The penultimate column of table III contains informa- 
tion about the purported superconducting gap (if appli- 
cable) observed by STM in the given experiment. In the 
samples for which bulk superconductivity is expected, 
STM detected at least one gap of the expected mag- 
nitude on most surfaces. An appropriate gap was de- 
tected on every surface with an ordered 2x1 struc- 
ture, most surfaces with y/2 x \/2, and even most sur- 
faces with disordered mixtures of these two structures. 
Massee explicitly observed "identical superconducting 
gaps on both y/2 x y/2 as well as 2 x 1 or other sur- 
face terminations" [TT4] . However, extreme disorder does 
seem detrimental to the observation of a superconduct- 
ing gap: Bover[90] observed no superconducting gap on 
very disordered surfaces which lacked atomic resolution, 
while Massee [TT4] observed no superconducting gap on 
very disordered surfaces with corrugation >2A. 

Therefore, according to the STM measurements of 6 
different groups[90j EH [Ml HHJ UTS 1333, Bal22 and 
Sri 22 surfaces of nominally superconducting composi- 



tion do display at least one superconducting gap of the 
expected magnitude, except in cases of extreme disor- 
der. The one puzzling exception is Zhang's failure to 
observe a superconducting gap on any of his supercon- 
ducting samples with y/2 x y/2 surface structure. It is 
also notable that none of the groups who have claimed 
As termination [56, 98] 1102] . have even looked at nomi- 
nally superconducting samples. Therefore, it is an open 
question as to whether an As-terminated surface, which is 
presumably polar and different in doping from the bulk, 
can support superconductivity. 

We now close the section on STM of 122 materials, 
and look ahead to the sections on 11 and 111 materi- 
als, which will turn out to be much more conducive to 
surface studies, with reliable, well understood, and likely 
bulk-representative cleaved surfaces. So we ask: is there 
anything more to be learned from STM on 122 materi- 
als? Given the challenges associated with their surfaces, 
in comparison to the 11 and 111 materials which follow, 
it may be tempting to label the 122 materials a waste of 
time, and move on. However, several new developments 
show that it is crucial that we understand the surfaces 
of the 122 materials, and when their properties can be 
trusted to be representative of the bulk. Two of the most 
intriguing new Fe-based superconductors are of the 122 
family: BaFe2(Asi_ iC P :c )2 and K IC Fe2- 2/ Se2. 

BaFe2(Asi_ :E P ;E )2 becomes a superconductor without 
charge doping: isovalent chemical pressure from P dis- 
rupts the magnetic order and allows the material to 
become superconducting with T c as high as 31K[86j. 
BaFe2(Asi_ ;E P :E )2 is therefore an ideal material in which 
to isolate and study the causes and effects of the mag- 
netic order, without the complications of changing carrier 
concentration, or disorder directly within the critical Fe 
plane. (This view has recently been questioned [HE].) 
Furthermore, BaFe2(Asi_ :E P :E )2 is the highest-T c P- 
containing Fe-based superconductor, with T c comparable 
to hole-doped 122 materials, and higher than electron- 
doped 122 materials. This is especially surprising given 
the evidence for sign-changing gap nodes in the OP of 
this compound |119[ I120| . Such nodes are expected to 
significantly lower the T^ |121] , 

Additionally, BaFe2(Asi_ :E P cc )2 may actually be useful 
to clarify the cleaved surface of the 122 compounds, be- 
cause the P atoms are known to be smaller and therefore 
lower (closer to the Fe plane) than the As atoms [T22] . If 
the cleaved surface is indeed the Asi-^P^ layer, the ex- 
pected concentration x of P atoms should be clearly vis- 
ible on the surface. In fact, Massee imaged the cleaved 
surface of BaFe2(Asi_ :E P :E )2 (x=0.32, T C =31K) but he 
did not see any signature of ~ 1/3 of the surface atoms 
sitting lower than the others, as would be expected for a 
Asi-^P^ top layer [97]. 

An even more intriguing newcomer is K0.78Fe1.7Se2, 
with T C =30K[123 , by far the highest in the Se subfamily 
of Fe-based superconductors. More generally, this family 
(K,Tl,Rb,Cs) £C Fe2- 2/ Se2 has been found to superconduct 
at ambient pressure with T c up to 4QK |124j . The ma- 
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Material 
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Cleave T 


Meas. T 


Term. 


Structure 


SC gap? 


Ref. 


parent Bal22 


PES 


300K 


300K 


Ba 


disordered 




de Jong[TT2l 


parent Bal22 


/V-LEED 


300K 


20K 


As 


3. 8 A tetragonal 




Nascimento 98 


parent Bal22 


STM 


20K 


4.3K 


As 


\/2 x \/2 


_ 


Nascimento 98 


parent Bal22 


STM 


80K 


5K 


As 


\/2 x \/2 


- 


LifTOO] 


parent Srl22 


STM 


300K 


40K 
& 200K 


Sr 


disordered 2x1 




Hsieh|92] 


parent Sri 22 


ARPES 


10K 


10K 
& 200K 


Sr 


2x1 


- 


Hsieh[92] 


parent Srl22 


STM 


warm 
& cold 


5K 


As 


y/2 X \/2, 2 x 1 


- 


Niestemski[T02] 


Ca(Feo.97Coo.o3)2As 2 


STM 


cold 


4.3K 


? 


2x1 


- 


Chuang 56 


( "parent" ) 
















K cc Sri_ x Fe 2 As 2 


STM 


10K 


5.3K 


Sr 


2x1, disordered 


yes 2x1, 


Bover[90] 


(T C =32K) 












no elsewhere 




Ba(Fe .9Coo.i)2As 2 


STM 


~ 25K 


6K 


Ba 


2x1 


yes 2x1 


Yin|9T] 


(T C =25K) 
















Ba(Fei_ cc Co x ) 2 As2 


STM 


300K 


4.2K 


Ba 


disordered 2x1 


yes everywhere 


Massee[TT3] 


0=0,0.07) 
















Ba(Fei_ x Co cc ) 2 As2 


LEED 


300K 


4.2K 


Ba 


disordered 2x1 




Massee[TT3] 


0=0,0.07) 
















Ba(Fe .93Coo.o7)2As 2 


STM 


< 80K 
& 300K 


4.2K 


Ba 


y/2 X y/2, 2 X 1, 
other 


- 


Massee[95] 


Ba(Fe .93Coo.o7)2As 2 


7V-LEED 


< 80K 
& 300K 


4.2K 


Ba 


2x1 


- 


Massee[95] 


Ba(Fei_ cc Co cc ) 2 As2 


STM 


80-100K 


4.5-24K 


Ba 


n/a 


yes except large 


Masseep"T4] 


0=0.04, 0.105) 












corrugation 




Ba(Fei_ cc Co x ) 2 As 2 


STM 


300K 


4.5-24K 


Ba 


n/a 


yes except large 


Masseep"T4] 


0=0.07) 












corrugation 




Ba(Fei_ cc Co x ) 2 As 2 


STM 


? 


? 


? 


\/2 x \/2 

— >• 2 x 1 as x t 


yes 


Pan[T03l 


Ba(Fei_ x Co cc ) 2 As 2 , 


STM 


120K 


5K 


Ba 


mostly "s/2 x \/2, 


yes 2x1, 


Zhang 94 


0=0-0.16) 










rare 2x1 


no elsewhere 




Ba(Fei_ x Co x ) 2 As 2 , 


7V-LEED 


<20K 


<20K 


Ba 


V^xV^ and 2x1 


_ 


van Heumen|104| 


0=0.05,0.085) 
















K x Sri_ cc Fe 2 As 2 


STM 


cold 


2.5K 


? 


no atomic res. 


yes, 2 gaps 


ShanpT5lfTT6l 


(T c = 38K) 
















Ba(Fei_ x Co cc ) 2 As 2 


STM 


300K 


6K 


? 


no atomic res. 


yes, 2 gaps 


TeaguepTf] 


0=0.06,0.12; 




in Ar 












T C =14,20K) 

















TABLE III. Summary of the surface structures reported for A122 compounds. 



terial is found to exist on the edge of an AF insulator 
phase which possibly originates from superlattice order- 
ing of the Fe vacancies |124j . Remarkably, ARPES mea- 
surements show that K0.8Fe1.7Se2 (T C =30K) is strongly 
correlated (with band renormalization by a factor of 2.5 
compared to LDA calculations) and lacks the hole-like 
FS at the T point which is present in other high-T c Fe- 
based superconductors |125j . This compound therefore 
rules out nesting as the sole scenario to explain Fe-based 
high-T c superconductivity. This ARPES result is credi- 
ble because the measured FS volume is found to satisfy 
Luttinger's theorem with the bulk carrier concentration; 
this is also very encouraging that STM studies of the sur- 
face this material will be fruitful windows into the bulk. 

Very recently, K ;E Fei_ 2y Se2 has been studied by STM 
by two groups. First, Li et at demonstrated MBE growth 
on a graphitized 6H-SiC(0001) substrate resulting in 



atomically flat (llO)-oriented films of K x Fei_ 2y Se2 [126 . 
STM imaging of these films demonstrated phase sepa- 
ration into stoichiometric superconducting KFe2Se2, and 
insulating K ;E Fei_ 2y Se2 in which the Fe vacancies tended 
to order into a y/b x y/b structure. In fact, a single iso- 
lated Fe vacancy in an otherwise stoichiometric region 
was shown to locally suppress superconductivity. 

Acomplementary STM study by Cai et al of the (001) 
face of cleaved bulk Ko.73Fei.67Se2 supported the idea of 
phase separation |127j . Cai et al found that the cleaved 
surfaces of multiple samples were predominantly unsta- 
ble and disordered due to a partial layer of mobile K. 
Occasionally, they found regions with an exposed, atom- 
ically ordered Se surface, and good superconducting spec- 
tra. However, despite extensive searching, they never 
found signatures of Fe vacancies in these regions, thus 
supporting Li's conclusion that the good superconduct- 
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ing regions are stoichiometric KFe2Se2- Furthermore, in 
these 'good' superconducting regions, Cai et al found a 
weak \/2 x y/2 reconstruction, suggesting coexisting an- 
tiferromagnetism. 

Both BaFe 2 (Asi_ a; P :c )2 and (K,Tl,Rb,Cs) a; Fe2- y Se2 
are likely to provide insights into the mechanisms behind 
Fe-based high T c superconductivity. As with any doped 
material, inhomogeneity must play a role, and therefore 
it is crucial to apply a local tool such as STM to resolve 
the mysteries of the DOS in these intriguing new com- 
pounds. We hope that insights from the cleaved surfaces 
of Co- and K-doped Ba,Sr,Ca-122 can be extended to al- 
low rapid characterization of surfaces, and deep insights 
about bulk DOS in these intriguing new materials. 



C. 11 Materials 

The simplest family of Fe-based superconductors, the 
'11' materials, have actually been studied for decades, 
but their superconductivity was discovered only after 
the sudden attention on the iron-pnictides [T28] . Be- 
cause there is no intervening ionic layer between the 
Fe(Te,Se,S) planes, there is no polar surface catastrophe, 
so the ll's were the first of the Fe-based superconductors 
which were ideally suited for STM study. 

Massee et al first showed that the surfaces of 
Fe1.07Teo.55Seo.45 appeared identical when cleaved at 
room temperature or low temperature, with every Se/Te 
atom visible, and no apparent surface reconstruction [95 . 
Cold-cleaved LEED (17K) experiments showed only the 
4A-periodic tetragonal spots with no extra reconstruc- 
tion. In this sample the 7% interstitial Fe were seen 
as very bright spots exactly between four visible Se/Te 
atoms, suggesting that they sat beneath the surface but 
contributed significant extra local DOS. 

Kato et al studied Fe1.05Teo.85Seo.15 with onset 
T C =14K[111 . They also saw the excess Fe as bright 
spots on the surface. Spectroscopy at T=4.2K showed 
A=2.3 meV and 2A/&bT c ~3.8, reasonable values which 
give no cause to suspect the surface is not represen- 
tative of the bulk. The standard deviation <ja=0.23 
meV gave <ta/A ~ 10%, in good agreement with a a/ A 
in Ba(Feo.gCoo.i)2As2 [91], and about half the relative 
gap variation reported in the cuprate superconductor Bi- 
2212 [129 . Kato et al did observe larger variations in the 
background DOS at higher energies, and speculated that 
this variation was due to the excess Fe. They further 
speculated that this inhomogeneity at higher energies 
may be responsible for the apparent ~10% inhomogene- 
ity in the superconducting A, but they did not address 
this point directly using spectra at the same locations 
above and below T c for normalization. 

Fridman et al reported atomic resolution STM on 
FeTei-^Se^, at nominal doping levels of x=0.3 and 0.5, 
with critical temperatures T c ^12K [T30] . The samples 
were warm cleaved in vacuum, then the spectra were mea- 
sured at 300 mK. Fridman confirmed the nanoscale inho- 



mogeneity in the topography and spectroscopy. Gener- 
ally, the spectra displayed a broad, V-shaped background 
out to ±10 mV, with a sharper gap structure A ~2-4 
meV, which vanished above T c . 

He et al reported STM on FeTei-^Se^ with x=0 and 
0.45, determined by energy dispersive x-ray spectroscopy 
(T C =14K for the £=0.45 sample) [I31j . The samples were 
cleaved and imaged at T=80K. They claimed a small sur- 
face corrugation of 8 pm (but this would be dependent on 
the resolution and work function of the tip) . They identi- 
fied the "bright" atoms as Te and the "dark" ones as Se, 
and noted that the heights were bimodal, separated by 
Az = 44.7 ± 12.0pm, and that the fraction of bright and 
dark atoms matched the known bulk Te/Se composition 
ratio. Both filled state and empty state topographic im- 
ages showed the same Az, supporting the conclusion that 
the apparent height difference was a real geometric effect, 
rather than a difference in DOS. This STM- measured Az 
was almost a factor of two larger than the value measured 
by x-ray scattering |132j . The larger value observed by 
STM could be a systematic experimental effect or a sur- 
face relaxation. He et al observed nanoscale clustering of 
the Se and Te, but noted that the spectra were homoge- 
neous (indistinguishable between Te and Se) out to ±100 
meV. They speculated that the length scale of chemical 
inhomogeneity was likely shorter than the superconduct- 
ing coherence length. They further speculated that the 
apparent chemical inhomogeneity of the surface Se/Te 
may have been balanced by the opposite species on the 
opposite side of the Fe plane. 

The previous experiments all studied the cleaved sur- 
faces of bulk Fei+^Tei-^Se^. Better control over the 
crystal composition can be obtained in films grown by 
molecular beam epitaxy (MBE) [50] . Song et al grew and 
performed STM imaging on FeSei +2 films, demonstrat- 
ing superconductivity for z < 2.5%. Sparse excess Se 
appeared as dimers and eventually, on increasing concen- 
tration, merged into an ordered y/h x a/5 surface structure 
with a 0.5 eV insulating gap. By optimizing the substrate 
temperature during growth, Song was able to reduce the 
number of defects to fewer than 1 in 70,000 atoms, and 
reliably achieve superconductivity. A study of supercon- 
ducting gap vs. film thickness showed that the T c scales 
inversely with the film thickness as T c (d) = T c o(l — d c /d) 
where T c q is the bulk T c of 9.3K, and d c is the mini- 
mum thickness for superconductivity, estimated as 7A, 
approximately two Se-Fe2-Se layers. 

These experiments addressed the structure and chem- 
istry of the surface of FeTei-^Se^. Two of them observed 
a superconducting gap of ~2-4 meV, but without suffi- 
cient energy resolution to make any concrete statements 
about pairing symmetry. Recently, two different studies 
by Hanaguri et al [27] and Song et al [57] addressed the su- 
perconducting pairing symmetry in the FeTei-^Se^ fam- 
ily, to be discussed in section (V) 
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D. Ill Materials 

The '111' family of Fe-based superconductors is also 
ideally suited for STM study. LiFeAs cleaves well with 
Li splitting evenly between cleaved surfaces. Calculations 
show no surface states in LiFeAs[133 . The first (unpub- 
lished) work on the 111 family shows an unreconstructed, 
atomically resolved surface, on which 6 different types of 
impurity states, as well as magnetic vortices have been 
imaged [49 . 

Very recent QPI imaging experiments on clean 
LiFeAs surfaces purport to demonstrated p-wave 
superconductivity [3 1, or to resolve the in-plane 
anisotropy of several coexisting gaps [134]. 



E. Other distant cousins 

Other Fe-based superconductors include Sr 2 V03FeAs 
(21311) with T c =37.2 K|T35j. (Sr 4 Sc 2 6 )(Fe 2 P 2 ) (42622) 
with T C =17K[136], and Sr 3 Sc2Fe 2 As205 (32225) with 
a possible T c ~20K |137j . In addition, there are 
several series with increasing Fe 2 As 2 interlayer sepa- 
ration: (Fe 2 As 2 )(Sr 4 (Sc,Ti) 3 8 ), (Fe 2 As 2 )(Ba 4 Sc 3 7 .5), 
and (Fe 2 As 2 )(Ba 3 Sc 2 5 ) with T c up to 28K[138]; also 
(Fe 2 As 2 )(Ca n+ i(Sc,Ti) n 2y with n=3,4,5 and with T c up 
to 42K for n=5 [139 . Many of these appear to have 
mirror planes which will be suitable for STM study, but 
so far they have not been successfully grown as large 
single crystals. For example, the largest so far in the 
Sr 2 V0 3 FeAs family are 300 fim [140] . smaller than is 
convenient for STM study. 



V. SUPERCONDUCTING ORDER 
PARAMETER 

A superconducting OP, which describes the symmetry 
of the pairing state, has both magnitude and phase. We 
start this section with an overview of theory and early 
experiments on the OP in Fe-based superconductors. We 
then explain the STM determination of the OP, starting 
with phase, which may be more relevant to the pairing 
mechanism, and proceeding to magnitude. 

Within the first few weeks of Hosono's discovery, ev- 
ery imaginable superconducting OP was proposed. Be- 
cause all five Fe d bands may cross the Fermi level, the 
number of possible OPs which can live on these five FSs 
is large. Some candidates were nicely summarized pic- 
torially by Hicks et a/ |141j . Because of the presence 
of Fe and suspected magnetic ground state of the par- 
ents, and the early LDA calculations showing nesting be- 
tween hole and electron FSs at the AF wavevector [58 , 
two early papers [142} [143 convincingly argued for an s=b 
OP (sketched in figure fl2^i). A spin- mediated pairing 
mechanism would require a sign flip between the phases 
of the paired carriers, forcing a plus sign on one nested 
FS and a minus sign on the other. In this scenario, there 



would be no need for a sign change on a single FS, i.e. 
no need for gap nodes. 

To cement a picture of spin-mediated s± pairing, three 
pieces of evidence are necessary: (1) demonstration of 
the absence of nodes on each FS; (2) demonstration of a 
phase flip between FSs; and (3) demonstration of a spin 
resonance at the nesting wavevector connecting the two 
opposite-signed FSs. 

Early ARPES studies did not show nodes in 
NdFeAsO . 9 F .i[69], Bao. 6 Ko.4Fe 2 As 2 [144] (see fig- 
ure 16 1), or Ba(Fei_ :E Co :c ) 2 As 2 [l45] (see figure [l6j b 
and cj. Many subsequent ARPES experiments have 
also failed to find evidence of nodes. But ARPES 
experiments do not provide phase information. 

Several phase-sensitive tunneling experiments have 
been performed. First, a search for half- flux-quantum 
vortices trapped in the native grain boundaries of poly- 
crystalline NdFeAsO0.94F0.06 (T C =48K) detected none, 
ruling out d-wave order in this compound |146j . How- 
ever, a detection of half-quantum magnetic flux in a 
loop between NdFeAsO0.8sF0.12 and conventional su- 
perconducting niobium [I47j was explained by frustrated 
Josephson coupling due to a sign change in the OP [148 . 
In combination, these two experiments strongly sug- 
gested s± symmetry in NdFeAsOi-^F^, but did not give 
proof for the pairing mechanism. A neutron scatter- 
ing experiment on another 1111 FeAs superconductor, 
LaFeAsOi-^F^, detected a resonant spin fluctuation that 
peaked sharply in the superconducting state, at the nest- 
ing wave vector [T49] . (The resonant spin fluctuation had 
already been seen in Ban ^Kn AFe9As9. |150j ) 

If we are willing to combine evidence from several com- 
pounds, a picture arises of a spin-mediated s=b OP in 
1111 and 122 materials. Nonetheless, other calculations 
have shown that various OPs with nodes on the FSs were 
not too different in energy, and may be obtained under 
some circumstances [T2T[ I151[ 152 j. Furthermore, there 
is a long list of experiments which showed evidence for 
nodal OPs in various Fe-based superconductors [T531I161] . 
Most of these experiments were sensitive only to the pres- 
ence or absence of low- lying quasiparticles, thus impuri- 
ties could give the same signature as a nodal OP. In the 
early days, as sample quality was still rapidly improv- 
ing, it was commonly believed that the elimination of 
impurities would eliminate the appearance of nodal su- 
perconductivity. 

However, the 11 family refused to fit the mold. In par- 
ticular, the parent FeTe compound did not show the same 
spin ordering at the nesting wavevector, but rather spin 
ordering at an alternative non-nesting wavevector [162 . 
Therefore, the idea of spin-mediated s± pairing in this 
material was called into question. As described in sec- 
tion [HDj STM is an ideal tool to take on this question. 
STM can simultaneously look for nodes (via the pres- 
ence of low-lying quasiparticles in dl/dV spectroscopy) 
and can perform a phase-sensitive test of the OP via 
QPI imaging. Several theoretical approaches predicted 
discriminating QPI patterns for expected OPs in the Fe- 
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based superconductors [25, 1631 1164j . 



A. Gap Symmetry 

s± order parameter in FeTe^^Se^^ Hanaguri et al 
performed STM experiments on FeTei-^Se^ with x ~ 
0.4 and bulk T c in the range 13K to 14.5K|27]. He first 
observed that the low temperature (T=400 mK) spectra 
were fully gapped, as shown in figure [13^. These fully- 
gapped spectra provide strong evidence against nodes in 
the OP. 

We address two possible arguments against the node- 
less gap: (1) STM sensitivity to a quasiparticle state 
ipk is exponentially suppressed with increasing in-plane 
momentum k[ll . One might therefore argue that STM 
lacks the sensitivity to detect quasiparticles arising from 
gap nodes far from the T point, e.g. on M-centered 
FSs. Two counter-arguments follow, (a) The FT-STS 
images (e.g. figure 12 i) show broad QPI peaks at q2 
and <?3, near the corners of the BZ, demonstrating the 
sensitivity of this STM and this particular tip to these 
large-/c quasiparticles. (b) STM spectroscopy at the 
same temperature, on a related material that is too 
clean to support impurity-induced quasiparticles, shows 
a V-shaped spectrum which could arise only from nodal 
quasiparticles [57 . These gap nodes are likely to live on 
the M-centered FSs (see figure p^), again suggesting 
that STM is generally sensitive to nodal quasiparticles 
even far from the T point. (2) Fridman et al found a 
V-shaped gap for FeTei-^Se^ crystals of a similar com- 
position, with nominal x ~ 0.3 (measured at at T=300 
mK)[130 . However, Fridman's samples were cleaved at 
room temperature, allowing the possibility of significant 
surface contamination. In fact, their images did not 
achieve the same clean atomic resolution as shown by 

It therefore seems likely that 



Hanaguri in figure 12 



Fridman's V-shaped gap arose as a consequence of scat- 
tering from surface contamination. 

Although Hanaguri's experiment provides strong evi- 
dence for a nodeless gap, it does not rule out the pos- 
sibility of gap anisotropy. Angle-resolved specific heat 
(ARSH) measurements on FeTeo.55Seo.45 suggest deep 
gap minima or nodes along the Fe-Fe direction (T-M di- 
rection in the BZ of figure comparison with 
Hanaguri's tunneling data, Zeng et al speculated that the 
gap minima or nodes detected by ARSH must live on the 
inner electron pocket (M pocket), i.e. at maximal k val- 
ues where the tunneling matrix element would be most 
severely suppressed, thus hiding the low energy quasipar- 
ticles from STM. However, Hanaguri's QPI data suggests 
that quasiparticles are highly visible to STM even when 
located at the M point, so reconciliation of ARSH with 
STM does not require that the nodes or gap minima live 
near the M point. Indeed, a closer look at Hanaguri's 
T = 0.4K spectrum shows energy difference ^0.75 meV 
between the lowest-energy quasiparticles, and the maxi- 
mum of the coherence peak. This difference greatly ex- 



ceeds the expected thermal broadening of ~ 4/c^T = 0.14 
meV, suggesting that the observed broadening is due in- 
stead to gap anisotropy. Therefore, Hanaguri's spectra 
support the likelihood of deep gap minima, but strongly 
argue against gap nodes. 

In the s± OP model, the electron and hole FSs have 
opposite signs, as sketched in figure [12^ l. The FSs are 
nested, with three possible nesting vectors sketched in fig- 
ure plfo. Two of these nesting vectors, qi and #3, are sign- 
preserving in the s± model, while the third, #2, is sign- 
changing. Upon application of a magnetic field, pinned 
vortices are expected to become a new source of sign- 
changing scattering [41 . Furthermore, magnetic impuri- 
ties, and non-magnetic resonant scatterers are both ex- 
pected to enhance sign-changing scattering and suppress 
sign-preserving scattering upon application of a magnetic 
field [42 . In the s± pairing model, all three mechanisms 
are therefore expected to enhance the coherence factors 
for qi and qs nesting, at the expense of ^2 nesting. This 



expectation is supported by the data shown in figure 12 i. 
The sign-changing s± OP is therefore strongly supported 
by STM in FeTe . 6 Se .4. 
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FIG. 12. (a) Schematic FSs of FeTeo.6Seo.4- Red square de- 
notes the unfolded BZ (one Fe per unit cell), while the dashed 
square at 45° denotes the tetragonal BZ. The hole FS at the 
r point (blue) and the electron FS at the M point (red) are 
thought to have different signs in the s± pairing scenario, (b) 
QPI will be dominated by the scattering vectors qi, ^2, and 
q3 as shown. Blue circles represent sign-preserving scatter- 
ing, whereas red diamonds represent sign-reversing scatter- 
ing, (c) Topography, showing clean atomic resolution (setup: 
V'sampie = — 20 mV; /set = 100 pA). Four bright spots corre- 
spond to extra Fe atoms, (d) Magnetic- field induced change 
in QPI intensities, Z(q,E 1.0meV,£ 10T) - Z(q,E 
l.OmeV, B — 0T), supports the s=b scenario. Figure from Ref. 

ET] 



A concern was raised by Mazin et a/ [166 , who pointed 
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out that all three q- vectors discussed by Hanaguri seem 
too sharp to be ascribed to QPI, and in fact each cor- 
responds to long range structural or magnetic order in 
the crystal. In particular (referring to the crystal lattice 
parameters defined in figures [7^ and b), q\ = 27r/aF e -Fe> 
q2 = 27r/ao, and qs = 2tt / 'ar would each correspond to 
long range crystalline order and give rise to sharp g-space 
peaks. In contrast, elastic scattering between each of the 



pairs of points on the FSs shown in figure [12ft would re- 
sult in broad q- space peaks, with diameters matching the 
diameters of the FSs, approximately 15-20% of the BZ. 
Mazin suggested that the observed intensity vs. magnetic 
field trends for q 2 and qs could instead be explained by a 
field-induced suppression of both the superconductivity 
and the spin density wave (SDW) and its assumed con- 
comitant surface reconstruction, which would enhance 
the structural Bragg peak at qs within the superconduct- 
ing gap energy and suppress the SDW-induced peak at 

Hanaguri countered 1 167] with linecuts through the q 2 
and qs peaks which showed that each peak was made 
up of two components with distinct energy and field 
dependence. Each peak could be separated into a central 
sharp Bragg peak, and a broader QPI peak with the 
expected width ~ 20% of the BZ. Hanaguri's s± gap 
symmetry conclusions were drawn only from the broader 
QPI peaks. Hanaguri argued that the coexistence of 
sharper Bragg peaks does not negate these conclusions. 

Nodal order parameter in FeSe There remain 
many credible claims of a nodal gap in Fe-based 
superconductors [T53Ul61| . However, from bulk studies, 
it can be hard to conclusively rule out impurity scatter- 
ing as the source of apparent nodal quasiparticles. Most 
recently, Song et al performed a beautiful set of exper- 
iments on the single purest Fe-based material studied 
to date: MBE- grown FeSe, with fewer than one defect 
per 70,000 Se atoms (i.e. no visible impurities in a 100 
nm square field of view) [57] . Clearly, there should be 
no impurity-induced low-lying quasiparticles in a sam- 
ple this clean. Nonetheless, at the lowest temperatures 
(down to T=0.4K), Song found a V-shaped gap (fig- 
ure [13)3) , which provided clear evidence of nodal super- 
conductivity. 

How could the same 11 material have two different 
OPs at different dopings? In fact, both cos k x cos k y 
and cos k x + cos k y pairing functions are consistent with 
the same point group symmetry, so they may naturally 
coexist [57]. The gap function is given by 

A(&) = Ai cos k x cos k y + A2(cos k x + cos k y ). (15) 



Figure 14 1 shows the nodal lines of these two terms in 
the unfolded BZ. Depending on the relative magnitudes 
of Ai and A 2 , the nodal lines of the full function A(fc) 
may or may not pass through the electron FSs (figure 
[14)3). It has been shown by functional renormalization 
group calculations that when the interlayer coupling t z is 
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FIG. 13. (a) dl/dV measured on FeTei-^Se* with x -0.4 
and bulk T C =14.5K (and apparent local T c ~11K). The low- 
est temperature dl/dV is clearly fully gapped, (b) dl/dV as 
a function of T for FeSe with T c — 8K. Although topographic 
imaging shows fewer than 1 in 70,000 defects, the low T dl/dV 
is still V-shaped, indicating the presence of low energy quasi- 
particles. In both Hanaguri's and Song's spectroscopy, the 
bias modulation amplitude was set to 0.1 mV rms . 



weaker, the nodal lines are more likely to intersect with 
the electron FSs, giving a d-wave OP[ 12Tj fl52| . 




/c x (1/a F e-Fe 



FIG. 14. Schematic BZ of FeSe, in the unfolded scheme, with 
one Fe per unit cell, (a) Hole- like FS around the F point, and 
electron-like FSs around the M points are shown in green. 
Nodes of the cos k x cos k y OP are shown in black. Nodes of 
the cos /c x +cos k y OP are shown in red. (b) Model showing the 
nodes of the combined Ai cos fee cosk y + A2(cos k x + cosk y ) 
OPs with Ai/A2=0.35. Increasing A1/A2 causes the blue 
line to move inwards, removing the intersections of the nodes 
with the M-centered electron FSs. 



In FeTei-^Se^, the Se atoms lie closer to the Fe plane 
than the Te atoms [T3T1 1132] , which may lead to reduced 
interlayer coupling t z . It is therefore consistent that the 
FeSe studied by Song would have lower interlayer cou- 
pling than the FeTeo.6Seo.4 studied by Hanaguri, and 
would have greater tendency to d-wave order. (It is also 
consistent with the finding that LaFePO, where P sits 
close to the Fe plane, has nodes [T68l I169| , whereas 1111 
materials with As, which sits farther from the Fe plane, 
seem not to have nodes [69 ). 

Therefore, low temperature STM on clean samples 
has proven to be a sensitive test of both the presence 
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of nodes, and the phase of the OP. It is important to 
note again that STM is presumed to be more sensitive 
to the near-r states [170] . Low-T spectroscopy on clean 
FeSe showed a V-shaped gap spectrum, although the 
gap nodes are thought to be on the M-centered FS. 
The fully-gapped spectra observed by STM at low T 
are therefore likely truly indicative of no nodes on any 
FS[27 . More generally, Song's work strongly suggests 
that STM is sensitive to the presence of nodes on any FS. 

p-wave order parameter in LiFeAs Recently, an- 
other QPI imaging study argued for p-wave supercon- 
ducting gap symmetry in LiFeAs [31]. LiFeAs with 
T C =18K was cleaved and imaged at low temperature, dis- 
playing an atomically resolved field of view with sparse 
impurities around which were centered clear interference 
patterns. The Fourier transforms of the interference 
patterns strongly resembled the band structure seen by 
ARPES [171], suggesting that the JDOS was dominated 
by a van Hove singularity at the T point. The authors 
computed the expected q- space QPI patterns for s±, 
p, and s + + OP symmetries. From a qualitative compar- 
ison, they suggested that p-wave symmetry was the best 
match, but they did not rule out more exotic symmetries 
such as s + id. However, their images showed some dis- 
tortions, and were not consistent with other STM reports 
on LiFeAs [491 [T34]. 



B. Gap magnitude 

The magnitude of the superconducting gap can be 
measured by bulk techniques such as heat capacity and 
electromagnetic absorption. The momentum-dependence 
of the gap can be directly measured by ARPES, which 
has so far been primarily limited to temperatures above 
4.2K (with rare exceptions [171]). ARPES measures only 
filled states, and cannot be used in a magnetic field. 
More precise energy resolution, as well as access to empty 
states and field-induced states, can be obtained through 
low T tunneling measurements [172]. If a clean surface 
can be obtained, STS can determine the gap with atomic 
spatial resolution, and with energy resolution limited 
only by the measurement temperature, which has been 
pushed down to 10s of milliKelvin |173j . STM energy res- 
olution is ~ 4/c#T, where the broadening stems from a 
smearing of DOS by ~ ksT on either side of the Fermi 
level, for both the sample and tip states. On some 
Fe-based superconductors, STS has been performed at 
temperatures down to 300 mK[l30] or 400 mK[27l 157]. 
Using QPI imaging [46 or clever fitting techniques j!74j , 
the /c-dependence of the gap magnitude can also be 
quantified by STM, although applications of these tech- 
niques to the Fe-based superconductors have not yet been 
published [134]. 

A typical superconducting dl/dV spectrum measured 
by STM has one or more pairs of roughly particle-hole 
symmetric peaks in the DOS, which may be interpreted 



as the superconducting coherence peaks. The simplest 
measure of the superconducting gap A is \ the energy 
between the two peaks. This technique has been used 
with success to analyze large dl/dV datasets to generate 
'gapmaps' which show the local variation of the OP mag- 
nitude in cuprate superconductors [175 and in Fe-based 
superconductors [9T| 1113] (see figure 17). 

In some cases, the coherence peaks may be fit very ac- 
curately to an expected functional form, originally pro- 
posed by Dynes et a/ |172j . Dynes noted that the super- 
conducting DOS p s (E) = \E\/VE 2 - A 2 could be gener- 
alized to take into account a finite quasiparticle lifetime 
by writing 



Ps(E,T) =Re 



E-iT 



y/(E - iT) 2 - A 2 



(16) 



where V is the inverse quasiparticle lifetime. Further gen- 
eralizing to an angle-dependent gap A(#), plugging into 
equation [4j and differentiating gives 



dl 
dV 



dE 



•Re 



de 
E-iT 



s=E-eV 



y/(E - irp - Mfff 



(17) 



However, in most cases, fits to equation 17 fail the 
common sense 'chi- by-eye' test (see figure 15). Possible 
causes for this failure include unaccounted-for energy- 
dependence of the background DOS, energy or momen- 
tum variation of the tunneling matrix element, or energy- 
dependent limitations to the quasiparticle lifetime [176 . 
The energies of the coherence peaks seem by eye to 
be reasonably well found by the various fits, but it 
must be noted that changes to the inverse quasiparti- 
cle lifetime V (also sometimes called the 'Dynes broad- 
ening parameter' [172 ) can dramatically shift the fit 
value of A. Additionally, thermal broadening causes 
the coherence peaks to appear to shift to larger energy, 
even as the gap A is actually closing upon increased 
temperature [177] . Therefore, in complex materials with 
multiple gaps and/or non-trivial momentum dependence, 
it is not clear whether careful fitting actually results in 
more accurate gap determination than the simple method 
of extracting half the distance between apparent coher- 
ence peak maxima (see figure 15 a and b). 

Fasano et al studied the 'as-grown' (air-exposed 
and uncleaved) surface of a ~100 jam single crystal 
of SmFeAsOi-^F^ with nominal x=0.2 and measured 
T C =45K[72 (figure [21^). Spectra, acquired from base 
temperature 4.2K, showed a V-shaped gap that disap- 
peared around the bulk T c . The half-distance between 
'peaks' (or kinks) in the spectra, gave an average gap 
value A p =7 meV, with standard deviation ~6% of the 
mean value. The reduced gap 2A P /ksTc ~ 3.6 was in 
good agreement with point contact spectroscopy on the 
same compound. However, after averaging (binned by 
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A p ), and normalization (division by a smooth polyno- 
mial fit to the high energy part of the spectrum above 



SmFeAsO-,, x F x (T C =45K) 
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the peaks), the spectra were fit to the Dynes equation 
with both s and d wave functional forms for A(#), yield- 
ing A s values from 4.2 to 4.8 meV (2A s /k B T c - 2.2), 
and Ad values from 5.1 to 5.8 meV (2 Ad/ ksT c ~ 3), as 
shown in figure [15^. Neither fit gave a reduced gap value 
close to the expected BCS values of 2A s //cbT c = 3.5[178 
or 2Ad/kBT c = 4.3 |179] , nor did either fitting method 
give consistent values of the inverse quasiparticle lifetime 
T. Because the gap values from <i-wave fitting were only 
~22% smaller than expected, in comparison to the s- 
wave fits which were ~36% smaller than expected, the 
authors concluded that their data supported d-wave pair- 
ing. 

Fasano et al also studied a 'dip- hump' feature in the 
spectra, which occurred just outside the superconduct- 
ing gap, and was interpreted as a signature of a collec- 
tive mode. The energy E 2 of this feature had an even 
larger spread around its mean value of ~ 15 to 20 meV. 
Unlike the inner superconducting gap, the outer dip- 
hump feature was not symmetric about the Fermi level. 
On the positive energy side (empty states), the energy 
of the presumed collective mode (Q = — A p ) was 
anti-correlated with the presumed pairing strength (A p ). 
Note that Q ranged from 2 to 8 meV in Fasano' s sam- 
ple, giving Q/ksTc ~ 0.5 — 2. This low energy of the 
purported 'spin resonance' Q in Fasano's work, may be 
cause for some concern, as the energy of the magnetic 
resonance has been shown to scale with T c in cuprates as 
Vt r - 5/cRT.[T80j. in Bal22 as Vt r = 4.3/c B T c [l50], and in 
LaFeAsOi-^F^ with near-optimal x = 0.082 and T c =29 
K as O r = 4.4/crX1~ |149] . In summary, all of the energy 
scales reported by Fasano (A a ±, A^, and Q) were alarm- 
ingly low compared to measured values from related ma- 
terials. It is possible that surface contamination played 
a role in these results. 

The expected reduced gaps 2A//cbT c for weak- 
coupling s-wave and d-wave BCS superconductors are 
3.5 and 4.3, respectively [178 1 HL 79] . In the cuprates, ex- 
periments which decoupled the pseudogap by coherent 
tunneling [181] or normalization [T82] have found reduced 
gap values between 6 and 10. Superconducting gaps mea- 
sured by STM in Fe-based materials are shown in ta- 
ble |IV| There is a wide distribution of reduced gap val- 
ues 2A//c#T c , which may arise from surface quality is- 
sues, unavoidable surface states even in clean surfaces, 
or method of determination of the gap, be it careful 
Dynes fitting or simple peak-to-peak measurement. In 
most cases, the reduced gap suggests that although Fe- 
based superconductors are in a strong coupling regime, 
they are not as strongly coupled as cuprates. 



C. Two Gaps 

Early ARPES measurements detected multiple FSs 
for the 122 family, and showed that these FSs may 
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FIG. 15. Dynes fits to single- and double-gapped spec- 
tra, (a) d-wave and s-wave Dynes fits to 3 representa- 
tive spectra from SmFeAsOi- x F x with nominal x—0.2 and 
measured T C =45K[72 . Both d-wave and s-wave curves are 
indistinguishable by eye, yet result in very different val- 
ues for the superconducting gap A. Neither are visually 
satisfying fits, and both result in reduced gaps 2A/kBT c 
more than 25% too low for a BCS superconductor. (b) 
s-wave fit to FeTen «sSen is pTT] also fails the 'chi-by-eye' 
test, (c) Fit using an extension of the Dynes form to two 
gaps, on Ba(Fei_ a ,Co x ) 2 As2 with x=0.06 and T.=14K[TT7]. 
(d) Fit using an extension of the Dynes form to two gaps, 
on Ba(Fei_ x Co x ) 2 As 2 with x=0.12 and T C =20K[TT3- (e) 
Fit using an extension of the Dynes form to two gaps, on 
Bai_ x K x Fe 2 As 2 with x=0.4 and r.=38K[TT5]. Without in- 
troducing a finite Dynes broadening parameter T, neither 
nodal nor nodeless gaps produce good fits, (f) Fit to the 
same sample in (e), after introducing a finite Dynes broad- 
ening parameter T. Still, neither nodal nor nodeless gap fits 
pass the common sense 'chi-by-eye' test. 



support two gaps of rather different magnitudes j!44j . 
For example, in optimally hole-doped Bao.6Ko.4Fe2As2 
(T C =37K), the two gaps are A=6meV and A=12meV, 
as shown in figure 16 1. However, in optimally electron- 
doped Ba(Fei.925Coo. 075)2 AS2 (T C =25.2K), these two 
gaps are much more closely spaced, at A=6.6meV and 
A=5meV |145] , as shown in figure [l6| b and c. 

The tunneling rate along the z direction is ex- 
pected to be strongly suppressed with increasing in- 
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TABLE IV. Superconducting gaps, measured by STM. 

Average gap values were determined by observing the coherence peaks; individual spectra were checked with s-wave Dynes 
fitting. 



plane momentum \k\\ |[TT 



Therefore, the tunneling 
matrix element of the M-centered bands is expected 
to be strongly suppressed with respect to that of the 
r-centered bands [T70] . In optimally electron-doped 
Ba(Fei.925Coo.o75)2As2, there is only one band centered 
at the T point, with A=6.6meV |145] , as shown in fig- 



ure [16j b and c. Therefore, it is not surprising that in 
extensive STM studies of atomically resolved surfaces of 
electron doped Ba(Fei_ :E Co :c )2As2[9TJ II 13] . the contri- 
bution from the outer bands has not been strong enough 
to show a clear second gap (with expected similar en- 
ergy A 7> 5 ~5meV) in addition to the observed gap of 
A ~6.5meV. 

However, two other very recent STM studies of non- 
atomically-resolved surfaces of superconducting Bal22 
have shown two gaps. First, Shan et al studied 
Bao.6Ko.4Fe2As2, the identical hole-doped compound pic- 
tured in figure [TBfc ,. In non-atomically-resolved regions, 
they found two gaps of A ~3.6meV and A ~8.1meV 
(fit values vary slightly with the choice of quasiparticle 
lifetime, as shown in figure 15 , c and d). 

Second, Teague et al resolved two gaps on the electron- 



doped Ba(Fei_ :E Co ;E )2As2, as shown in figure 15 , e and f. 
This study differed from all previous STM studies on the 
same material in that the sample was not cleaved in vac- 
uum (rather in Ar at room T), and no atomic resolution 
was obtained, even after cooling to low T. This gives a 
clue that the imaged surface may be contaminated, so the 
tunneling electrons may have their momenta scrambled 
by passage through the contaminant layer. This might 
allow enough mixing between bands to circumvent the 
rule of thumb that the tunneling matrix element for the 
M-centered bands will be suppressed [1 70 . 



Finally, a recent follow-up study on the same 
Bao.6Ko.4Fe2As2 sample[115 showed evidence of three 
gaps at ^3.3 meV, ~7.6 meV, and ~10 meV, on the 
non- atomically resolved surface [TT6] . In conclusion, two 
or more gaps have been resolved only in STM measure- 
ments which lack atomic resolution, suggesting that sur- 
face disorder scrambles the quasiparticle momentum. 
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FIG. 16. Fermi sheets and superconducting gaps of Bal22. 
(a) Schematic of Fermi sheets and superconducting gaps 
of hole-doped Bao.6K .4Fe2As2 (T C =37K), as measured by 
ARPES fTH] , (b) Schematic of Fermi sheets of electron- 
doped Ba(Fei.925Coo. 075)2 As2 (T C =25.5K), as measured by 
ARPES [145] . (c) Isotropic superconducting gap values of the 
Fermi sheets shown in (b). 




FIG. 17. (a) Gap map recorded by Yin et al over 

a 20x20 nm 2 area of Ba(Fe .9Co .i)2As 2 at 6.25 K[9T]: 
and (b) by Massee et al over a 18.9x18.9 nm 2 area of 
Ba(Feo.93Coo.o7)2As2 at 4.2 K |113j . In both cases, the re- 
ported gap values are computed as half the distance between 
coherence peaks. 



Here we discuss the spatial and energetic signatures of 
electronic ordered states of the parent compound. One 
spatial signature, detected by STM, is a rotational sym- 
metry breaking from C4 — >• (72. A much-debated en- 
ergy signature in the parent compound is a depression in 
the DOS near the Fermi level, which has been framed in 
the language of a possible 'pseudogap', in analogy to the 
cuprates. 



D. Gap Inhomogeneity 

A key advantage of STM over other gap measure- 
ment techniques is access to the gap variation on the 
nanoscale. In most Fe-based superconductors, modest 
nanoscale variations have been found in the gap magni- 
tude, with standard deviations amounting to ^10% of the 
average value A, as listed in table IV and exemplified in 
figure [17} In cuprates, superconducting gap variances of 
<ja/A >20% have typically been rep or ted [129]. More re- 
cently, it was shown that these gap variations in cuprates 
are heavily influenced by the pseudogap |182j . but even 
after the pseudogap is removed by normalization or in- 
dependent fitting, the superconducting gap variations in 
cuprates appear larger than those in the more weakly 

A corre- 



coupled Fe-based superconductors [182] 1185] 
lation between the superconducting gap variance oa/A 
and the reduced gap 2A/&sT c has been suggested in Fe- 
based superconductors [97 . 



VI. PARENT COMPOUND 

Understanding the mechanism of superconductivity re- 
quires an understanding of the 'normal state' of the 
parent compound out of which superconductivity arises. 



A. C4 — >> C2 symmetry breaking 

Several theories predict C4 — » C2 electronic symmetry 
breaking in the parent compounds of Fe-based supercon- 
ductors, due to orbital ordering[186-189 or a nematic 
state [152} 11901 1191] . A nematic state is a liquid which 
breaks long-range rotational symmetry without breaking 
long-range translational symmetry [1 92 . Therefore, if the 
rotational symmetry is broken first by the orthorhombic 
crystal structure, the detection of C2 electronic order 
does not strictly constitute a nematic state. However, if 
the crystal symmetry is broken by only ~1%, as is typi- 
cal in the orthorhombic phase of these materials [193], an 
electronic C4 — >• C2 symmetry breaking of 50 — 200% or 
more has been argued to be too large to be solely due to 
a pre-existing structural effect, and has been ascribed in- 
stead to an independent electronic nematic state [5 6] [57]. 

Some theories have predicted incipient magnetic order- 
ing above the crystallographic orthorhombic transition at 
temperature T s |190[ I191j (see figure [6^ for a representa- 
tive phase diagram). According to these theories, above 
T s , there are only fluctuations of the sublattice magnetic 
order, which will eventually, on cooling, lead first to the 
structural transition and then the full magnetic transi- 
tion. Below T/v, the AF magnetic order additionally 
breaks translational symmetry, so the phase is no longer 
strictly a nematic. Between T s and T/v, the C4 symmetry 



23 



2n/a Fe . 




R ao detwinned (A) /g\ 

R t twinned (A) ^ 

R ao detwinned (B) 

R t twinned (B) 


R b0 part, detwinned (C) (f) - 

' R, twinned (C)/| >^ 


/ ? 0-75 
/ | 0.70 




// I"' 9 






160 170 180 


150 160 170 180 



0.008 

of 



50 100 150 200 250 50 100 150 200 250 

Temperature [K] 

FIG. 18. STM, ARPES, and transport experiments have all 
provided evidence for C2 symmetry in 122 materials, (a) QPI 
images on 3% Co-doped Cal22[56]. The dominant wavevec- 
tor corresponds to length ~ 8ciFe-Fe (shown in red). From 
inspection of twin boundaries, the authors determine this to 
be along the orthorhombic ao axis (the AF axis), (b) ARPES 
momentum space image (second derivative) from 3.5% Co- 
doped Cal22, at OmeV binding energy [194]. There is a nest- 
ing vector between parallel FS sections at q = 27r/8.3aFe-Fe 
(shown in green), (c) QPI dispersion from Chuang 56 . (d) 
Transport anisotropy measured on detwinned Bal22 [T95] . 
The high conductivity direction in Bal22 is the orthorhom- 
bic a direction (the AF axis), (e-f) Transport anisotropy 
measured on detwinned Cal22 [T96] . (e) The resistance de- 
creases along the strain direction (with dominant a-oriented 
domains), (f) The resistance increases along the orthogo- 
nal direction (with dominant 6-oriented domains). Taken to- 
gether, these results confirm that in Cal22 as well, the high 
conductivity direction is the orthorhombic a direction (the 
AF axis). 



is broken but translational symmetry is not yet broken, 
fitting the definition of a true nematic phase [190 . 

Because as-grown crystals are composed of multiple 
orthogonally-oriented orthorhombic domains, intrinsic 
anisotropy is challenging to measure via bulk experi- 
ments. Local STM experiments were the first to show 
evidence for large C4 — >• C2 electronic state symmetry 
breaking [56]. Chuang et al reported QPI measurements 



from cleaved Ca(Feo.97Coo.o3)2As2, at the underdoped 
edge of the superconducting dome (shown as a red ★ in 
the phase diagram of the similar Ba(Fei_ ;E Co :E )2As2 sys- 
tem in figure^). As expected at that doping, no signa- 
ture of superconductivity was seen in STM spectroscopy 
at T=4.3 K. In these non-superconducting parent com- 
pounds, deep in the orthorhombic phase, Chuang high- 
lighted the radically different electronic behavior along 
the two orthorhombic axes (recall the orthorhombic unit 
cell in figure ^f)p) . 

Along the a-axis (the longer, AF axis), the DOS 
showed self-similarity on a length scale d ~ 8aF e -Fe- The 
authors emphasized the existence of a 'static' electronic 
nano-object of size d, which folded the band structure by 
a wavevector q a = 27r/rf |197] , labeled in figure 18 1. 

Along the 6-axis (the shorter, ferromagnetic axis), 
the authors noted significant QPI dispersion, labeled as 
wavevector q^ in figure [TH^i. The dispersion of q^ along 6, 
in contrast with the 'static' wavevector q a ~ 27r/8ape-Fe 
along a, led the authors to posit a nematic band, which 
'evolves along one axis only'. 

Chuang showed excellent quantitative agreement be- 
tween the q\y dispersion along 6, and the dispersion of 
strongly nested sections of the T-centered ot2 hole band, 
observed by ARPES [T98] in undoped Cal22 (a rigid band 
chemical potential shift was used to account for the dif- 
ferent dopings in the STM and ARPES experiments). 
However, the ARPES measurement were carried out on 
twinned Cal22, giving the appearance of C4-symmetric 
nesting. 

Given the C2-symmetric surface structures already 
demonstrated on various 122 materials [90 1 [9TJ |95j 1102] , 



it is important to ask whether the STM-observed C2- 
symmetric QPI results from a surface or bulk phe- 
nomenon. Chuang did a meticulous job to separate the 
C4 — )> C2 symmetry breaking of the observed QPI from 
the more obvious 2x1 reconstruction of the surface, 
which also breaks C4 symmetry, but at a 45° angle (see 
topography in figure 10 z). To detect the very small 
1% length distortion and 1° angle distortion of the or- 
thorhombic phase, the authors scanned the same field 
of view in two different directions, and found both sur- 
face reconstruction twin boundaries (see red arrows in 
figure [lOfc), and orthorhombic twin boundaries (see green 
arrows in figure [Tofc), at distinct locations. The authors 
therefore empirically linked the QPI anisotropy to the 
bulk orthorhombic orientation, rather than the surface 
reconstruction orientation. It is clear from this work 
that (a) STM can image the bulk electronic structure of 
122 materials, and (b) bulk electronic structure is signif- 
icantly more anisotropic than one might naively expect 
from the 1% lattice orthorhombicity. 

Although the experimental data gave clear evidence for 
bulk C2 symmetry, theoretical objections were raised to 
the claim of a nematic origin for this C2 symmetry [60 . 
Mazin objected that by definition, nematic symmetry 
doesn't apply at Chuang's low temperature where the 
long range AF order breaks translation symmetry. Fur- 
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thermore, Mazin insisted that a nematic band is not 
needed in order to explain Chuang's observations. Us- 
ing a few empirical values to simply parameterize a 
plausible SDW-reconstructed FS, Knolle et al calculated 
the QPI patterns due to scattering from pointlike non- 
magnetic impurities and magnetic impurities oriented 
along the z axis [199]. Knolle's calculations reproduced a 
C2-symmetric QPI pattern with wavevectors which were 
quantitatively similar to Chuang's measured q a and 
but unlike Chuang's data, the calculated QPI patterns 
dispersed along both orthogonal axes. Mazin et al used 
a different approach, starting from ab initio FS (shown 
in figure [7^) [60 . Mazin emphasized that no quantitative 
comparison should be possible without accurate knowl- 
edge of the surface doping and of the scatterers. He there- 
fore made no assumptions about the scattering potential 
and computed only the JDOS rather than the full QPI 
pattern. His calculated JDOS showed C2 symmetry, em- 
phasizing a dispersive scattering vector along the b axis, 
but lacking any features, dispersive or not, along the a 
axis. Therefore, both Knolle and Mazin were able to 
reproduce some of the qualitative features of Chuang's 
measured QPI, in particular the C2 symmetry, without 
invoking a nematic band [60) 1199] . 

Following Chuang's local STM measurements, a break- 
through by Chu et al allowed majority detwinning of the 
bulk orthorhombicity in 122 materials (for example, ~5 
MPa uniaxial pressure resulted in 86% volume fraction 
of one orientation) |195j . Chu reported that detwinned 
Ba(Fei_ cc Co ;E )2As2 developed in-plane electronic trans- 
port anisotropy near the structural transition. The re- 
sistivity along the shorter, ferromagnetic axis, p^, be- 
came greater than p a . The resistive anisotropy reached 
a maximum value of ~2 for compositions near the un- 
derdoped edge of the superconducting dome (see fig- 
ure [l8|i). Even for temperatures well above the or- 
thorhombic structural transition, uniaxial stress induced 
a resistivity anisotropy, suggesting a substantial nematic 
susceptibility. 

Tanatar et al verified the transport anisotropy in 
detwinned Bal22, and extended the measurement to 
detwinned Cal22, allowing direct comparison with 
Chuang's STM data |196j . Tanatar showed that upon de- 
twinning, resistivity decreased along the orthorhombic a 
axis but increased along the orthorhombic b axis in both 
compounds. In both materials Tanatar found that the 
resistive anisotropy was largest at T s with Pb/p a ~ 1-2 
in Cal22 and ~1.5 in Bal22. For Cal22 this anisotropy 
was observed only below T s , suggesting a first order tran- 
sition, and diminished upon further cooling, reaching 
about 1.05 at T ~50 K and remaining constant at lower 
temperatures. For Bal22 the anisotropy was observed 
both below and above T s , suggesting a second order tran- 
sition with nematic fluctuations above T s . 

Following Chu's demonstration of bulk detwinning, 
ARPES experiments on detwinned samples showed the 
anisotropic band structure of Cal22 [T94] and Bal22[200 . 
In Cal22, Wang et al confirmed the presence of a strongly 



nested FS, with nesting vector approximately match- 
ing Chuang's QPI vector (figure plfo). Furthermore, 
the detwinned ARPES augmented Kondo's work 1 198] by 
showing that the nesting was (72-symmetric. In Bal22, 
the FS showed elongated regions of the same orientation 
and location, but without the strong nesting of Cal22, 
so a susceptibility to nematicity was emphasized only at 
temperatures above T s [200 . 

Some aspects of these experimental results from STM, 
transport, and ARPES seem counterintuitive, or even 
contradictory. First, one expects that conductivity along 
the ferromagnetic direction (b) should be higher than 
along the AF direction (a). But transport experiments 
consistently show that p^ > p n |195[ 1196] . Second, an 
a-oriented nematic band as posited by Chuang to ex- 
plain the observed anisotropic QPI would be expected 
to lead to p a > pb in contrast to the observed trans- 
port anisotropy. Third, given the dramatic anisotropy of 
the band structure, both from the experimentally mea- 
sured strongly nested FSs |194] , and theoretically cal- 
culated fully broken symmetry between the Fe xz and 
yz orbit als [196], one might expect a dramatic electronic 
transport anisotropy. In fact, only weak to moderate 
anisotropy up to a factor of 2 is observed. Finally, the 
maximal electronic transport anisotropy in Cal22 was 
found right below the structural transition T s rather than 
at lower temperatures where the obvious C4 symmetry 
breaking OPs (either magnetic moment or structural or- 
thorhombicity) would be expected to increase [196 . 

To qualitatively resolve some of these issues, Mazin 
noted that, while the calculated FSs completely break 
the C4 symmetry, which is fully reflected in Chuang's 
QPI images, the individual pockets are very three dimen- 
sional, so that the calculated conductivity is comparable 
for all three directions [60 , and the anisotropy may take 
either sign according to parameter details to be worked 
out [196]. 

Most recently, follow-up doping-dependent STM stud- 
ies by Allan et al [197 have suggested a more complete 
and quantitative explanation for the counterintuitive and 
apparently contradictory results of STM, transport, and 
ARPES experiments. Allan et al showed that in the pres- 
ence of C2-symmetric electronic structure, Co dopants 
either establish or pin static electronic dimers, oriented 
along the a axis, with a length scale of d ~ 8aF e -Fe- He 
showed that these dimers may act as extended scatter- 
ing objects which, in conjunction with the nested band 
structure measured by Wang et a/ |194] , lead to QPI 
which is quantitatively consistent with that observed 
by Chuang, and to resistive anisotropy consistent with 
observations [195, 196 . 

Several other parent compounds of the Fe-based super- 
conductors have subsequently been studied by STM. So 
far, all have shown evidence for a C2 ordered state. 

Li et q/ |100] claimed evidence of C4 — ^ C2 symmetry 
breaking in Bal22 from subtle differences between the 
shapes of 45° and —45° domain walls, as shown in figure 
[TO^ i. However, unlike Chuang et al, Li did not image a 
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FIG. 19. C2-symmetric QPI in LaFeAsO. (a-b) 20 nm 
square topographic images from the two inequivalent cleaved 
surfaces of LaFeAsO. (a) FeAs surface (setup: Sample = 
100 mV; /set = 50 pA). (b) LaO surface (setup: Sample = 
— 100 mV; /set = 40 pA). (c) Fourier transform of a 40 nm 
square dl/dV map of the LaO surface at -75 meV. The gray 
circles at the corners show the Bragg peaks of the La surface 
atoms, 27r/aLa-La where OL a -La ~4A. (d) Dispersions of the 
QPI signatures at q\ and q^. Figures taken from Ref. 1751 



twin boundary in the bulk orthorhombic lattice, and thus 
could not demonstrate that the subtle domain wall orien- 
tation differences were coincident with the orthorhombic 
domain. Therefore, it remains possible that this subtle 
observation may be due to surface state anisotropy or 
even to tip anisotropy. 

STM images of the LaFeAsO parent compound also 
showed C2-symmetric QPI on one of the two inequiv- 
alent cleaved surfaces, shown in figure 19, a and b[75]. 
There are a number of significant differences in the ac- 
quisition and interpretation of the QPI data from Cal22 
and LaFeAsO. The surface showing QPI in LaFeAsO was 
believed to be LaO (based on an observed increase in 
imaged surface impurity states on doping from parent 
LaFeAsO to LaFeAsOo.gF .i), while the imaged surface 
in Cal22 was thought with less certainty to be FeAs |1 10] . 
The imaged QPI in LaFeAsO was believed to originate 
from a surface state, based on theoretical predictions [79] 
and on the fact that the atomic topography was almost 
completely obscured. The imaged QPI in Cal22 was 
shown to originate from the bulk, based on the fact that 
it changes orientation coincidently with a twin boundary 
in the bulk orthorhombic lattice but not the 2x1 sur- 
face reconstruction. It should be noted that in LaFeAsO, 
unlike in Cal22, there is no apparent reconstruction of 
surface atoms themselves into a C2 symmetric state. 



Zhou observed QPI at long and short wavevectors, 
but chose to focus on the more-easily-quantifiable long 
wavevectors q\ and q2 (figure [l9fc), because they are 
easier to distinguish from the long-wavelength inhomo- 
geneity that results simply from the random impurity 
distribution. Chuang removed the effects of this long- 
wavelength inhomogeneity by subtracting a Gaussian 
centered at the g-space origin. Chuang focused on the 
shorter wavevectors, which are labeled q a and q^ in fig- 
ure us*. 

Zhou found qi and q^ to disperse over a large en- 
ergy range from -200 to +200 meV, in perfect synchrony, 
as shown in figure [l9] i. In contrast, Chuang reported 
that qb dispersed strongly, while q a did not disperse at 
all, as shown in figure [TBfc . Although Zhou's q\ and q<i 
wavevectors dispersed identically, the shape of their q- 
space weight distribution broke C4 symmetry. In Zhou's 
data, q\ appeared as an isolated spot, whereas q2 ap- 
peared elongated across almost the entire BZ. Chuang's 
data also showed QPI peaks elongated along the b axis. 
The main difference in the conclusions of the two papers 
stems from Chuang's claim that q a did not change with 
energy, in contrast to Zhou's changing q\ and <?2 along 
both orthorhombic directions. This claim of static q a 
led Chuang to posit a nematic band, in a melted stripes 
picture [201 . However, anisotropic QPI may arise simply 
from a scattering potential with C2 symmetric structure 
factors [30], as shown empirically by Allan et a/[197 . 

Further investigations of C4 — >> C2 symmetry break- 
ing were carried out on MBE-grown FeSe films with 
near-perfect Se surfaces (fewer than one Se vacancy per 
^70,000 sites) [57] . Song et al found that magnetic vor- 
tices were elongated by a factor of ~2 along one Fe-Fe 
direction with respect to the other. Furthermore, they 
imaged a twin boundary which separated two regions 
of orthogonally-elongated vortices, thus proving that the 
C4 symmetry breaking was not a tip effect. Song also 
saw faint stripes in the background DOS, with wave- 
length ~ 16ap e -Fe along the same direction as the vortex 
elongation. These stripes are reminiscent of the 8ap e -Fe 
stripes observed in Cal22[56 , but with approximately 
twice the wavelength. Although Song did not investigate 
the energy dependence of these faint stripes, he specu- 
lated that they arose from intra-pocket nesting of the T- 
centered hole pocket, which may be distorted by orbital 
ordering [T86HT881 [202] . 

Song also observed single atom impurities whose spa- 
tial resonances locally break C4 symmetry on two differ- 
ent length scales. The short length scale 2-fold symmet- 
ric shapes ('atomic-scale dimers'), similar to those seen 
in LaFeAsO [75] (figure [19^), may be explained simply 
by the ~50% probability occupation of each of the two 
geometrically inequivalent Fe sites in the unit cell, with- 
out requiring a deeper mechanism of electronic symme- 
try breaking. Indeed, Song et al show the atomic-scale 
dimers to be randomly oriented on both sides of their 
orthorhombic twin boundary. However, a longer length 
scale electronic dimer structure also surrounds the same 
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impurities; this ~ 8nm dimer is always oriented con- 
sistently with respect to the orthorhombic crystal direc- 
tion. In comparison with the orientation of the ~ 8ap e -Fe 
electronic dimers demonstrated by Allan et a/ [197] , it 
seems reasonable to identify the vortex elongation direc- 
tion (longer as the crystalline a-axis (the AF axis). 



B. Pseudogap 

Although the Fe-based superconductors are exciting 
materials in their own right, it must be admitted that a 
significant component of the interest is due to their role 
as a foil for the higher- T c cuprates. One of the most puz- 
zling unsolved problems in cuprates is the origin of the 
'pseudogap' - a suppression in the DOS near the Fermi 
level at temperatures far exceeding the superconducting 
Tn |2Q3j . It is therefore a natural question: is there a pseu- 
dogap in the Fe-based superconductors? The answer is 
still controversial, with pseudogap sightings claimed by 
NMR[204], ARPES[205], femtosecond spectroscopy [206], 
and transport experiments [207]. In parent Srl22, the 
SDW transition occurs at Tsdw ~190K, and opens a 
gap of ~60 meV detected by optical spectroscopy [208 . 

The first question in evaluating these claims is, what 
do we mean by a 'pseudogap'? In the cuprates, the 
term arose to describe a mysterious gap of unknown ori- 
gin, which appeared at temperatures far above the ex- 
pected superconducting gap. If we therefore take 'pseu- 
dogap' to mean a 'gap of unknown origin', then STM 
has measured many pseudogaps in the Fe-based super- 
conductors. If we restrict ourselves to discussing repeat- 
able gaps of unknown origin, then there does seem to 
be a frequently observed phenomena, manifesting as a 
sometimes-symmetric set of kinks around ± 20-50 meV, 
in a background V-shaped DOS. However, there is not 
yet evidence from STM that this phenomenon turns on at 
any specific temperature higher than T c . The few papers 
which did track the probable pseudogap energy range 
through the superconducting T c showed no evidence of a 
gap persisting above T c [TT4J [183] . 

The situation in Fe-based superconductors is in some 
sense the inverse of the situation in cuprates. In cuprates, 
there was a clear gap opening at temperatures far above 
T c , but no obvious electronic ordered state, supercon- 
ducting or otherwise, to take the blame for that gap. 
This led to over a decade of speculation about electron 
correlation effects, numerous hidden OPs, and the pos- 
sibility of phase incoherent precursor pairing far above 
the bulk superconducting transition[180, 203 . However, 
in Fe-based superconductors, there are two clear phase 
transitions following a line in phase space which is sug- 
gestively similar to the T* pseudogap in line cuprates. 
First, there is a structural transition from tetragonal to 
orthorhombic at T s . Second (or in some cases simulta- 
neously) there is a magnetic transition to a SDW state 
(typically an itinerant collinear antiferromagnet) at T/y. 
This has launched a search for the 'pseudogap' which 



should be associated with these transitions. This seems 
a bit of a misnomer, because the phase transitions are 
known, so if their gap is indeed detected, there will be 
nothing 'pseudo' about it. 

In addition to the known SDW phase transition, there 
is an even more basic reason to expect a depressed DOS 
over a broad energy range close to the Fermi level. The 
Fe-based superconductors are compensated semimetals, 
which means that the Fermi level is located just below 
the top edge of one band and just above the bottom edge 
of another, meaning that the DOS at the Fermi level itself 
will be low compared to the DOS farther into either of 
these bands. So we should expect that all DOS measure- 
ments will be superimposed on a roughly V-shaped or U- 
shaped background, which will have a minimum roughly 
at the Fermi level for a compensated semimetal. All high 
energy DOS curves have shown this background, as ex- 



emplified in figure 20 



In contrast, a gap which opens up due to a new elec- 
tronic ordered state (such as superconductivity or charge 
density wave (CDW) or SDW) will conserve states, push- 
ing states out towards the gap edge, causing the appear- 
ance of gap edge peaks or coherence peaks. So the 'pseu- 
dogap' we seek should be more than just a depression 
in the DOS; it should open up at a particular T and 
it should also show some pile-up of states at the gap 
edge. Because it is superimposed on the semimetallic 
background DOS, the pile- up may manifest as a shoul- 
der rather than a peak. 

Several studies have shown evidence for a 'pile-up' in 
the DOS, above the superconducting energy scale. This 
higher-energy pileup has been interpreted as a second 
superconducting gap (figure [l5j c -f)[lT5] 1117] , the sig- 
nature of a collective mode (figure |21p ,) [72 , or a surface 
state (figure [21)3) [102] . In some cases, a higher energy 
pile-up of states has been observed, but not interpreted 
(for example, see figure 13^) [27]. 

Some authors have explicitly labeled their higher en- 
ergy features as 'pseudogaps'. In Fe1.05Teo.85Seo.15, Kato 
et al saw a broad, asymmetric V-shaped DOS out to at 
least ±500 meV, with a kink around +300 meV, which 



he identified as a 'pseudogap feature' (figure 20i) [lll| . 
In lightly Co-doped Cal22, Chuang showed a V-shaped 
'pseudogap' spectrum out to ±100 meV, with a metal- 
lic (non-zero) DOS at Ep (figure |2lfc) [56 . In K-doped 
Bal22 with T C =37K, Wray showed a 'pseudogap' of ~50 
meV (figure |2l]i) [184 . Other representative examples of 
DOS features at 'pseudogap'-like energies can be seen in 



figure 2l3-h. 

Despite a suggestive title, "Pseudogap-less high 
T c superconductivity in BaCo x Fe2- a; As2" , Massee did 
not definitively rule out a pseudogap in Co-doped 
Bal22 [TT4] . What he did demonstrate was that, unlike in 
the cuprates [182] , the apparent ~10% spatial variation in 
the energy of the superconducting gap cannot be blamed 
on the presence of a pseudogap. Using the same normal- 
ization technique as Boyer et a/ |182j , Massee divided the 
spectra below T c by the spectra at the same location just 
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FIG. 20. High energy dl/dV curves measured on several 
Fe-based superconductors. All curves are consistent with a 
roughly compensated semimetallic DOS, displaying a broad 
miminum near the Fermi level, (a) dl/dV from a series of 6 
compounds Ba(Fei- cc Co cc )2As2, where x is the nominal con- 
centration (from starting element ratio in the crystal growth 
process) [94] . All curves show a pronounced feature at ap- 
proximately -200 meV, which has been explained as a surface 
state, from calculations and ARPES measurements 104j. (b) 
A higher resolution spectrum from Ba(Fei_ cc Co cc )2As2 with 
nominal x=0.1 and T C —25K shows the low-energy supercon- 
ducting gap as well 105 . (c) Spatially averaged DOS of the 
parent Sri 22 compound 102 . This average was taken from 
a region displaying the \[2 x \[2 structure, but is qualita- 
tively similar to spectra from regions displaying the 2x1 
structure, (d) Spatially averaged DOS, measured at T=4.2K, 
from FeTeo.85Seo.15Hn]. 



above T c . Unlike Boyer, Massee found that the spatial 
variation in the apparent superconducting gap (defined 
as the distance between coherence peaks) was not re- 
moved by normalization. However, all 3 samples studied, 
with #=0.08, 0.14, and 0.21, ranging from underdoped to 
optimally doped, did show a relative depression of DOS 
near in the non-superconducting state spectra. This 
depression was most pronounced in the most underdoped 
sample, as shown in figure [21} . 

The spin-density-wave transition in 

Ba(Fei_ x Co :c )2As2 ranges from Tjy ~135K in the 
parent compound, down to ~20K where the phase 
boundary enters the superconducting dome and actually 
folds back slightly at doping x = 0.06 (see figure [6^) [53] . 
If the purported 'pseudogap' in this system were a 




FIG. 21. Spectra showing pseudogap-like features, (a) 
Spatially averaged DOS as a function of temperature from 
SmFeAsOi-xFcc with nominal x = 0.2 and T C =45K [72 . Low 
energy feature is identified as the superconducting gap; higher 
energy feature is identified as a collective mode. (The phase 
diagram of SmFeAsOi-xF^ shows no overlap between the 
SDW and superconducting states, as shown in figure^ from 
Ref. I5ip . (b) A linecut along the parent Srl22 shows a ubiqui- 
tous feature at around ±15 meV |102j . (c) An average spec- 
trum from Ca(Fei_ cc Co cc )2As2 with measured x — 0.03 56 . 
This sample is expected to be in the orthorhombic, non- 
superconducting SDW state, (d) Hole doped Bai_ x K cc Fe2As2 
with T C =37K shows a superconducting gap marked by black 
arrows, and a second higher energy gap- like feature marked 
by blue arrows [TBI], (e) BaFe2As2 shows asymmetric fea- 
tures of unknown origin [98]. (f-h) Ba(Fei_ x Co x )2As2 from 
Ref. 94 (f) Parent BaFe2As2 in the orthorhombic state. 

(g) Ba(Fei_ cc Co cc )2As2 with nominal x = 0.1 and measured 
T C =25K. The measured T c indicates that this sample is not 
in the orthorhombic SDW phase (see Fig ^ and Ref. [53}. 

(h) Ba(Fei_ x Co cc )2As2 with nominal x = 0.16 and measured 
T C =9K, unambiguously not in the orthorhombic SDW phase. 

(i) Underdoped Ba(Fei_ x Co cc )2As2 with measured x = 0.04 
and T C =14K [114] . This sample is expected to be in the 
orthorhombic SDW phase with Tsdw ^70K. Spectrum ac- 
quired at T=20K shows the expected disappearance of the 
superconducting gap, but no evidence of a 'pseudogap' within 
the expected mean- field energy range for temperature Tsdw- 



mean-field gap due to the SDW, then one would expect 
the pseudogap to range from ~20 meV in the parent 
compound, down to ~3meV where it cuts off at x = 0.06. 
In particular, the expected pseudogap would be ~10 
meV in Ba(Feo.96Coo.o4)2As2 studied by Massee, but 
he saw no evidence of a pseudogap edge out to ±50 
meV[TT4]. 
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In the end, what matters about any DOS feature iden- 
tified as a 'pseudogap', is its relationship to supercon- 
ductivity. If it is the signature of a competing state, or a 
collaborating state, then it is worth understanding how 
it interacts with superconductivity in momentum space. 
However, in these materials, particularly the 1111 and 
122 materials which do not cleave nicely, there are likely 
numerous surface states which may also manifest as fea- 
tures in the DOS. In our enthusiasm to draw connections 
with the cuprates, we must be cautious not to confuse 
these features with a pseudogap. 



VII. VORTEX STATE 

Magnetic vortices are important both for their rel- 
evance to technological applications (vortices must be 
well-pinned in order to allow high critical current J c ), 
and for their scientific interest as nanoscale windows into 
the low T normal state. The spatial shape and energetic 
signature of vortex core states may be used to deter- 
mine the pairing symmetry. The size of vortex core states 
may be used to measure the superconducting coherence 
length £. STM has now imaged vortices in 122[9T | [TT5] , 
in ll[27j|57], and in 111 compounds [49]. 



A. Vortex Pinning 

One of the very exciting early findings in Fe-based su- 
perconductors was the unexpectedly high critical field 
Hr-o |2Q9j and strong native vortex pinning[210 . A sum- 
mary of the technologically relevant properties of Fe- 
based superconductors, two years after their discovery, 
can be found in a review by Putti |211j . 

Although critical current is a bulk property, measur- 
able by bulk techniques, the quest to improve critical 
current relies on the pinning of vortices, typically by the 
judicious introduction of nanoscale defects which locally 
depress superconductivity. The system may then find the 
lowest energy state in which vortex cores (which locally 
destroy superconductivity over the length scale of the co- 
herence length £) are co-located with the defects, in order 
to avoid paying twice the energetic cost-of-destruction 
of superconductivity. The mechanism of vortex pinning 
may be investigated in aggregate by bulk critical current 
measurements, but in the case of unknown distributions 
of pinning sites, a local probe may be used to investigate 
the pinning sites of individual vortices [212 . 

The first STM images of vortices in an Fe-based super- 
conductor, optimally electron-doped Ba(Feo.gCoo.i)2As2 
(figure [22^i) [91] showed a disordered array, indicating 
that vortex pinning forces trumped the inter- vortex in- 
teractions which would drive vortices to form an ordered 
lattice. In the presence of strong pinning sites, the vor- 
tex arrangement with respect to those sites depends on 
the anisotropy of the material. In highly anisotropic su- 
perconductors, a one-dimensional vortex line may split 



like a stack of pancakes into point-like objects with the 
freedom to move independently in each superconduct- 
ing layer [2 13 . In this scenario, pancakes may find pin- 
ning sites independently in each layer, resulting in a high 
correlation between observed vortex and pinning loca- 
tions in any given layer. This scenario is realized in 
Bi2212[214 . In a more isotropic superconductor with 
strong pinning sites, the vortices must remain as line ob- 
jects, and can bend only slightly between layers to max- 
imize their overlap with point impurities throughout the 
bulk. In this scenario, there may be very little observ- 
able correlation between vortex locations and impurities 
in any single layer. This latter scenario was observed by 
Yin in Ba(Feo.gCoo.i)2As2 (figure 22i)[3 , giving evidence 
for strong pinning and low anisotropy in electron-doped 
Bal22. In contrast, the vortices in hole-doped Bal22 
were seen to form a hexagonal lattice, indicating weaker 
pinning in this material (figure 22 3) |115j . This supports 
the possibility that the Co atoms themselves, doped di- 
rectly into the FeAs layer of Bal22 (unlike the K atoms 
which are doped into the Ba layer), may act as the strong 
pinning sites in electron-doped Bal22. 

As would be expected, in clean FeSe (fewer than 1 
defect in 70,000 Se sites), vortices were also seen to form 
a hexagonal lattice (figure [22^ ) [57] . 

Unpublished conductance images by Hanaguri et al 
on nominally stoichiometric LiFeAs show a disordered 
vortex arrangement, indicating strong pinning [49 . How- 
ever, the topographic images of Hanaguri's LiFeAs sur- 
face showed a native impurity concentration correspond- 
ing to ^ 2% of the Fe sites. 

The superconducting coherence length £ is related to 
the upper critical field by H C 2 = <£o/27r£ 2 . In optimal 
electron-doped Bal22 (T C =25K), Yin et al measured a 
coherence length of £ = 2.9nm, corresponding to an 
upper critical field of H c2 = ^ /2tt£ 2 = 43T[91]. In 
optimal hole-doped Bal22 (T C =38K), Shan et al mea- 
sured a coherence length of £ = 2.2nm, corresponding to 
H C 2 = 75T |115j . Therefore, electron-doped Bal22 has 
stronger pinning at moderate fields, whereas hole-doped 
Bal22 seems to have a larger upper critical field. 



B. Vortex Core States 

In a conventional s- wave superconductor, theory pre- 
dicts and experiment confirms that the destruction of su- 
perconductivity in the vortex core will result in quasipar- 
ticle bound states at energy h A 2 /g^ |215j , which appear 
as a zero bias peak in the DOS [216 that splits into two 
symmetric peaks and eventually merges into the coher- 
ence peaks on moving away from the vortex center [2T7] . 
In d-wave cuprate superconductors, particle-hole sym- 
metric subgap states have been observed with energies 
approximately ± A/4 [2131 I2T8] . 

Yin et al found that in electron-doped Bal22, the 
vortex core destroyed the superconducting gap and 
coherence peaks, leaving behind a larger V-shaped 
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FIG. 22. (a) 106 nm square image of dl/dV at the ap- 
proximate coherence peak energy, -5 meV, showing vortices 
in Ba(Fei_ cc Co cc )2As2 with nominal x — 0.1 and T C =25K, as 
imaged at H=9T and T=6K (setup: V^ampie — —5 mV; 7 se t — 
10pA)[91 . Vortices are the broad blue regions, depressions 
in the DOS at the coherence peak energy. Near-surface im- 
purities are also visible as sharper, darker depressions in the 
DOS at the coherence peak energy. Vortices show no ordered 
lattice, and no correlation with the locations of surface im- 
purities, (b) 130 nm x 50 nm image of zero bias conduc- 
tance, showing vortices in Bao.6K .4Fe2As2 with T C =38K, as 
imaged at H=9T and T—2K (setup: V^ampie — 100 mV; 7 se t — 
200 pA) [115| . In this case the vortices show a hexagonal lat- 
tice with slight disorder . (c) 60 nm square image of zero bias 
conductance, showing vortices in FeSe with T c ^8K, as im- 
aged at if=8T and T=0.4K (setup: Sample 10mV;/ S et 
100 pA) [57] . Here vortices show a hexagonal lattice with slight 
disorder. Notably, they also show marked anisotropy, at least 
a factor of 2 x difference in size along the orthorhombic a and 
b axes, while the lattice distortion itself is only 0.5%. 



background [91]. But the cores lacked any of the sub-gap 
peaks predicted or observed in conventional s-wave or d- 
wave superconductors. This experimental null result was 
verified by S. H. Pan, but has not been published (cited 
as a private communication in Ref lll5j) . 

In contrast, STM studies of optimally hole-doped 
Bal22 showed a near-e^ peak in the DOS at the vortex 
center, which split and merged with the coherence peaks 
away from the vortex [TT5] (figure 23 i). However, the vor- 
tex center peak was slightly offset from the Fermi level, 
as the authors suggest might be expected in the quantum 
limit in which T/T c < 1/ (/cf£) [219| . This phenomenol- 
ogy is very similar to that seen long ago in NbSeo fZLT] . 
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FIG. 23. (a) dl/dV spectra along a 104A trajectory through 
the center of a single vortex in Bao.6Ko.4Fe2As2, at i7=9T, 
and T=2K [TT5] . (b) and (c) dl/dV spectra, spaced 2nm apart 
on a-axis and 6-axis trajectories through the center of a single 
vortex in FeSe at if=lT and T=0.4K[57]. 



In unpublished work on LiFeAs [49 , Hanaguri also 
showed a vortex core peak in the DOS, offset from sf, 
which split and evolved towards the coherence peaks 
away from the vortex center. Presumably, LiFeAs too 
was in the quantum limit at the 1.5K measurement tem- 
perature. Furthermore, the vortices in LiFeAs appeared 
star-shaped, with arms extending along either the Fe-Fe 
bond or diagonal direction, depending on the energy. 

STM studies of FeSe showed a vortex core state exactly 
at the Fermi level. Along one orthorhombic axis, the peak 
split, moved to higher energy, and eventually merged into 
the coherence peaks on a linecut away from the vortex 
center, as shown in figure [23|d. Behavior along the other 
orthorhombic axis was less straightforward, shown in fig- 
ure [23J3. In FeSe, the coherence length £ appeared much 
so the measurements were probably not in the 



longer, 

quantum limit where the vortex core state would be off- 
set from zero energy. 

A number of theoretical efforts have aimed to predict 
and explain vortex core states in Fe-based SCs, and par- 
ticularly to discern whether the observed structure of the 
core states can be used to elucidate the pairing symme- 
try. 

Vortex state calculations were performed on a 2-band 



model in the d T 



and 5^.2^2 pairing states 



and 



the Bogoliubov-deGennes equations were solved self- 
consistently by exact diagonalization. Resonance core 
states near the Fermi energy were found for both d and s 
pairing states. For the d x i_ y i pairing, the states are lo- 
calized while for the 5^2 ^2 pairing symmetry, they can 
evolve from sharp, energy-asymmetric localized states 
broad into extended ones with varying electron filling fac- 
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tor. To explain the non-observation of core states by Yin 
et al in Ba(Fei_ :E Co ;E )2As2, Hu et al show that the res- 
onant state is suppressed by an SDW, which may be a 
global coexisting phase (as in the underdoped materials) 
or may be locally induced by the vortex itself. 

Similarly, Jiang et al showed by self-consistently solv- 
ing the Bogoliubov-de Gennes equations based on a two- 
orbital model, that the magnetic field can induce a SDW 
state [221 in and around the vortex cores. Within the 
more favorable s± pairing, Jiang et al found that there 
are 2 kinds of vortices, those which induce SDW order, 
and those which don't, separated by a critical value of the 
Hund's coupling J c . When present, the SDW state ex- 
pands the vortex cores slightly, and suppresses the near- 
er core state. 

Araujo et al further considered the structure of a 
single vortex in two different 2-band models with s± 
pairing[222 . In a continuum model, they find a vor- 
tex core state which peaks at the Fermi level, while in 
a tight-binding model they find that the peak deviates 
from the Fermi level by an energy which depends on the 
band filling. Interestingly, they find that an impurity 
located outside the vortex core has little effect on the 
LDOS peak, but an attractive impurity close to the vor- 
tex core can almost suppress the vortex core state and 
shift the peak to negative energy. 

Wang et al calculated the vortex-core electronic struc- 
ture with both in-phase s-wave (s++) and antiphase s- 
wave (s+_) pairing, using four different tight-binding 
models[223 . They found a bound state LDOS peak at 
the core center, which splits away from the center and 
eventually merges to the gap edges, in qualitative agree- 
ment with experimental data from Ban aKn AFe9As9 |115| . 
They found that the sign of the peak energy uo p is pos- 
itive/negative if the Fermi level is near the bottom/top 
of the electron/hole bands, but that the peak energy is 
insensitive to the relative phase of the pairing gaps. The 
observed bound state in Ban ^Kn AFe9As9 pT5] is there- 
fore compatible, although not exclusively, with pair- 
ing. According to their calculations, the bound state 
should also be observed in Ba(Fei_ :E Co ;E )2As2, but they 
attributed its non-observation [9Tj to the Co dopants di- 
rectly in the FeAs layers, which may act as in-plane scat- 
tering centers, smearing out the bound state. (Later 
authors [224] noted that Wang's models were not self- 
consistent, and that the chemical potentials used for 
all 4 tight-binding models corresponded to undoped or 
electron-doped compounds.) 

Gao et al investigated the vortex states in 
Bai_ x K :c Fe2As2 by solving Bogoliubovde Gennes 
equations based on a phenomenological 2-band model 
with extended s-wave superconductivity and a com- 
peting SDW state 221 . Their result for the optimally 
doped compound without induced SDW is in qualitative 
agreement with the STM data on Ban ^Kn AFeoAso lTTS] . 
Similar to Hu et a/ |220| and Jiang et al [221] , they find 
that the SDW (present in underdoped samples) will 
reduce the intensity of the in-gap peak in the local DOS 



and transfer the spectral weight to form additional peaks 
outside the gap. Furthermore, the SDW will slightly 
enlarge the vortex cores, and reduce the rotational 
symmetry of the states from C4 to C2. Like Jiang et al, 
they speculate that an SDW state may be induced in the 
vortex cores of Ba(Fei_ :E Co ;E )2As2, thereby suppressing 
the vortex core states. 

In an effort to intepret Hanaguri's images of 
anisotropic vortex core states in LiFeAs[49 , the shape of 
the vortex core states was calculated by Wang et a/ [225 . 
One expects the vortex core states to decay exponentially 
with length £o = vf/ttA, where vp is the Fermi veloc- 
ity and A is the gap amplitude. In the direction of a 
gap node, the decay will follow a power law. However, 
Wang found that the geometry of the vortex core state is 
strongly influenced by the normal state band structure. 
If the gap itself is not highly anisotropic, the anisotropy 
of vf dominates, preventing direct observation of super- 
conducting gap features. 



VIII. CONCLUSIONS & FUTURE 

STM has made significant contributions to the under- 
standing of Fe-based superconductors. Spectroscopy and 
QPI imaging have provided a phase-sensitive determina- 
tion of the superconducting OP. After several years of de- 
bate about the existence of nodes in the superconducting 
gap, STM has provided clarification by demonstrating 
the evolution of the OP from nodal to nodeless upon dop- 
ing a single material [2 7 1 57[ . STM has also shown a C4 — >• 
C2 symmetry breaking in the parent and superconduct- 
ing states, highlighting the role of alternative and possi- 
bly competing electronic orders in these materials [56) 175]. 
Finally, STM has directly imaged the vortex state, as a 
route to better understanding and applying these excit- 
ing new superconductors^^ 57j [9TJ I115[ 1116] . Following 
these accomplishments, there remain several important 
open questions, to which STM is an ideal tool to provide 
answers. 

What is the nature of the C2 symmetry and its re- 
lation to superconductivity? QPI imaging and analysis 
throughout the entire BZ should be coupled with careful 
determination of the local orthorhombic orientation [56J. 
Most Fe-based superconductors are non-stoichiometric, 
and therefore may have inhomogeneity in their electronic 
anisotropy and strength of the superconducting state. 
Song's elongated vortex cores give a local measure of co- 
herence length anisotropy [5 7 . Real-space images of elec- 
tronic 'nematogens', and QPI-derived information about 
their scattering properties may give another measure of 
local anisotropy |56j I197j . STM should be used to mea- 
sure and correlate both the degree of local electronic 
anisotropy and the local superconducting energy gap. 

How does superconductivity interact with spin order? 
Spin-polarized STM (SP-STM) should be used to image 
local structures such as vortices or impurities which may 
pin or disrupt the spin order. Measures of broken (74 
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symmetry should be combined with spin-polarized scan- 
ning tunneling spectroscopy, to understand the energy 
scales of this broken symmetry state. 

How are vortices pinned in these materials? Larger 
area studies will be needed, with better identification of 
local impurities. 

We must also continue local studies of some 
of the most tantalizing new materials such as 
(K,Tl,Rb,Cs) iC Fe2- 2 ySe2, which challenge three-year-old 
beliefs of Fe-based superconductivity, such as the ne- 
cessity of nesting, and the deleterious effects of strong 
magnetic moments in the superconducting state. Bulk 
investigations of this material have been confusing due 
to possible inhomogeneity and difficulty controlling the 
ordering of the non-stoichiometric Fe. Conflicting re- 
ports claim bulk coexistence of superconductivity with 
antiferromagnetism[226], or phase separation |124j . 

In conclusion, Fe-based superconductors are an ex- 
traordinarily rich group of materials. Compared to the 
cuprates, their metallic ground state, relative malleabil- 
ity and electronic isotropy (typically less than a factor of 
two), large upper critical field and strong native vortex 
pinning may all make them more useful for applications 



in bulk |2 11] or thin film form [227]. Three years after the 
discovery of T C =26K in LaFeAsOi-^F^, new materials 
are still being discovered at a rapid rate [228 , which may 
lead to higher T c yet in Fe-based materials. They are in- 
triguingly similar to the higher-T c cuprates, and serve as 
a foil which has hastened the discovery of broader under- 
lying principles across even farther flung families of un- 
conventional superconductors [229 . With multiple pair- 
ing symmetries within the same subfamily, as well as mul- 
tiple temperature scales for magnetic sublattice ordering, 
structural transition, AF ordering, and superconductiv- 
ity within the same material, Fe-based superconductors 
are a complex playground in which to understand how 
various electronic and spin orders may support or inhibit 
superconductivity. 
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